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ABSTRACT 


This  work  represents  the  culmination  of  efforts  supported  by  Ballistic 
Research  Laboratories  and  Mobility  Equipment  Research  and  Development  Command 
to  develops  computer  simulation  system  for  large  scale  terrestrial  scenes. 

The  previous  efforts  have  produced  independent  models  for  solid  targets , cam- 
ouflage nets,  terrain  and  vegetation  ground  cover.  The  present  effort  has 
produced  automatic  generation  of  vegetation  input,  improvements  in  the  effi- 
ciency of  the  vegetation  model  and  increased  flexibility  of  the  terrain  model 
and  has  welded  the  independent  models  into  a single  system  capable  of  producing 
shadowed  photographs  of  all  or  any  part  of  a composite  scene.  The  system  has 
been  successfully  desionstrated  by  the  production  of  a high  quality  photograph 
of  a simulated  test  site. 
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1.  INTRODUCTION  AND  OVERVIEW  OF  THE  SYSTEM 

1.1  Background 

The  computer  system  described  in  this  report  represents  the  successful 
culmination  of  efforts  of  Mathematical  Applications  Group,  Incorporated,  sup- 
ported by  Ballistic  Research  Laboratories  and  Mobility  Equipment  Research 
and  Development  Cotmnand  and  others,  to  develop  a practical  tool  for  realistic 
computer  simulation  of  large  scale  terrestrial  scenes. 

The  previous  efforts  have  produced  independent  models  for  solid  targets*, 
camouflage  nets2,  terrain3  and  vegetation  ground  cover3 Although  independent, 
these  models  were  mutually  consistent  in  that  they  produced  images  in  the  same 
format  (block  format  image)  and  could,  therefore,  be  merged  into  a composite 
scene.  Associated  with  these  models  were  a group  of  utilities  for  merging,  mag- 
nifying and  shifting  of  images  and  a conversion  program  to  produce  from  the  block 
format  the  string  format  required  for  the  final  phases  of  picture  making. 

From  a practical  point  of  view,  there  were  a number  of  deficiencies  in  some 
of  the  models  and  in  the  system  as  a whole.  Specifically! 

a)  The  ray  tracing  elements  of  the  vegetation  model  required  too  much 
computer  time.  This  was,  in  part,  a result  of  the  multi-stage  ray  tracing  and, 
in  part,  a result  of  the  great  complexity  required  for  close-up  photography,  but 
wasted  on  distant  trees.  Thus,  both  a more  efficient  ray  tracing  was  suggested 
and  a different,  less  complex,  model  was  needed  for  the  low  resolution  situation. 

b)  The  generation  of  "forest"  data,  i.e.,  locations,  sizes  and  orientations 
of  the  trees  of  a scene,  was  performed  by  hand  and  this  placed  a literally  im- 
possible burden  on  the  user. 

c)  The  terrain  model  was  a modification  of  the  camouflage  net  which  tied 
the  color  and  texture  boundaries  to  the  resident  net  patterns.  This  gave  little 
flexibility  in  the  description  of  terrain.  Furthermore,  the  input  was  precisely 
the  input  for  a camouflage  net,  and  thus  did  not  proceed  directly  from  the  field 
data  available  to  modelers  of  terrain. 
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d)  The  generation  of  a composite  photograph  required  many  separate 
problem  set-ups  and  submissions:  one  for  each  modal  and  for  its*  shadow  and, 
for  the  vegetation  model,  separate  runs  for  each  tree  type.  The  calendar  time 
for  generation  of  such  a photograph  might  be  many  days. 

All  of  these  problems  have  been  addressed  successfully  and  the  solutions 
are  described  in  detail  in  subsequent  sections  of  this  report.  In  this  section, 
we  will  outline  the  improvements  and  guide  the  reader  to  the  more  detailed 
descriptions. 

1.2  More  Efficient  Ray  Tracing  for  the  Vegetation  Model 

A feature  of  both  the  construction  and  the  ray  tracing  of  the  vegetation 
model  was  the  multi-stage,  prototype  concept.  In  essence,  this  means  that  the 
tree  is  built  up  in  stages,  using  copies  of  a prototype  secondary  branch  to  con- 
struct a primary  branch  prototype,  and  copies  of  the  primary  branch  prototype 
to  construct  the  framework  of  the  tree.  In  ray  tracing,  the  process  is  reversed: 
when  a volume  containing  a copy  ie  struck  by  a ray,  the  ray  is  transformed,  i.e. , 
translated  and  reoriented,  into  the  (possibly)  magnified  world  of  the  prototype 
for  detailed  ray  tracing  within  that  structure.  Because  the  volumes  occupied  by 
the  copies  inevitably  overlap  a good  deal,  this  procedure  tends  to  invoke  much 
more  rry  tracing  than  is  required  for  the  more  open  structure  of  a one  stage 
model. 

Thus,  the  essential  feature  of  the  new  models  is  th*  elimination  of  some 
of  the  multi-staging  in  the  ray  tracing.  The  construction  of  the  tree  proceeds 
as  before  (from  the  same  input  data),  but  only  the  digitized  leaf  and  twig 
structure  is  represented  oy  copies  of  BOXes  in  the  system  of  the  tree.  All 
• woody  structure  - secondary  branch  base,  primary  branch  framework,  and  stem  - 

is  represented  by  truncated  right  cones  in  the  tree  system.  A three  dimensional 
lattice  covering  the  containing  box  of  the  tree  permits  analytic,  ordered  se- 
lection of  bodies  for  investigation  along  the  ray.  This  reduction  of  the  role 
of  snjlti-staging  has  resulted  in  a computar  time  saving  of  a factor  of  four 


or  five 


The  two  new  models  that  implement  the  now  ray  tracing  are  referred  to  as 
the  high  resolution  model  and  the  low  resolution  model.  The  high  resolution 
model  is  geometrically  identical  with  the  old  model.  The  low  resolution  model 
diffeis  from  the. high  resolution  model  in  that  the  BGXes  containing  leaf  and 
twig  structure  are  probed  "statistically",  i.e. , penetration  of  a struck  BOX  is 
probabilistic,  as  is  the  decision  whether  leaf  or  twig  has  been  struck.  Thus, 
this  model  differs  from  the  others  both  geometrically  and  in  its  ray  tracing 
method.  Because  another  "stage"  has  been  removed,  the  ray  tracing  is  roughly 
twice  as  fast  as  fcr  the  high  resolution  model.  The  statistical  data  for  this 
model  is  derived  from  ray  tracing  "experiments"  on  the  high  resolution  model. 

These  two  models  were  implemented  by  the  government  and  plied  to  four 
tree  types.  The  resulting  eight  data  sets  constitute  the  "tree  library"  de- 
scribed in  Section  2.  The  details  of  the  construction  technique  are  given  in 
Appendix  A. 

1.3  The  Data  Base 

The  data  required  to  describe  any  given  test  site  consists  of  five  libraries. 
Two  of  these  are  related  to  vegetation  ground  cover  and  are  1)  the  collection 
of  high  and  low  resolution  models  for  each  of  the  tree  types  of  the  region  and 
2)  the  forest  library.  Other  libraries  are  the  terrain  library,  the  Combina- 
torial Geometry  target  library  and  the  camouflage  net  library. 

Since  a prime  objective  of  this  effort  was  to  reduce  the  user  burden,  it 
was  decided  that  these  libraries  should  exist  on  a single  logical  unit  in  a 
highly  processed  state  and  be  accessed  by  "knowledgeable"  retrieval  routines, 
i.e.,  routines  which  are  aware  of  scene  limits  without  unnecessary  user  prompt- 
ing. The  degree  of  processing  in  the  stored  data  was  determined  for  each  model 
individually.  For  new  models  (vegetatic;.  and  terrain)  the  stored  data  are  the 
fully  processed  contents  of  arrays  and  common  blocks,  ready  for  ray  tracing, 
while  fcr  the  forest  description,  the  twelve  parameters  of  the  containing  BOX 
for  each  tree  are  stored.  The  net  and  CG  target  data  are,  however,  stored  as 


P— * 111  1 mp» H'.uwi1 1 1 1 '■!»»  till  ig nioy»u» ' u nw-ii»  i iimm  » ”i' ' "r  ■»•  iw™'-  - -» 

required  by  these  working  models:  c^rd  images  for  the  CG  and  a mix  of  card 
images  and  unformatted  restart  data  for  the  nets. 

iormat  of  this  tape  and  the  several  retrieval  methods  is  discussed 
in  detail  in  Section  2.  Data  selection  from  the  forest  library  and  data  or- 
ganizing for  the  vegetation  ray  tracing  will  also  be  discussed  there. 

1 , 1 The  Forest  Generator 

A stand-alone  program  for  the  generation  of  the  forest  library  has  been 
written  and  utilized  in  this  effort.  This  program  uses  resident  distribution 
functions  with  use  -specified  means  and  variances  to  distribute  trees  of  given 
types  within  given  boundaries  or  along  given  contours.  Amonc  the  properties 
selected  from  distribution  functions  are  location,  orientation,  height  and 
cluster  si~e. 

This  program  is  fully  described  in  Section  3,  which  also  contains  a user's 
guide  a.id  input  description. 

1.5  The  Terrain  Model 

Modifications  of  the  terrain  model  have  proceeded  by  separating  the  "in- 
stallation" procedure  from  the  ray  tracing  procedure,  and  building  a stand- 
alone program  to  produce  the  processed  data  required  for  the  data  base.  This 
stand-alone  program  requires  elevation  data,  arbitrary  color  boundaries  in  a 
horizontal  plane  and  associated  region  number  assignments.  Road  data,  as  in 
the  old  model,  are  also  required.  Only  minor  modifications  were  required  in 
the  ray-tracing  to  identify  the  color  boundary  within  which  a strike  occurs. 

The  »***tn llation  program  and  the  modifications  in  the  ray  tracing  are  fully 
described  in  Section  4.  That  section  also  contains  a user's  guide  ard  input 
♦ description  for  the  installation  program. 
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1.6  The  Image  Generating  System 

The  ’mage  generating  elements  and  shadow  generating  elements  of  the  four 
tasic  geometry  models  have  been  integrated  into  a system  which,  operating  from 
the  data  ba.“.;,  •-•an  produce  all  or  any  part  of  a scene  as  a block  format  image 
Ifor  furt  .ar  processing)  or  as  a string  format  image,  ready  for  the  color  pass 
of  the  M/.'-X  SYNTHAVISION  process.  The  other  components  of  the  system  are  the 
"scene  se  .ting"  components,  for  reading  and  processing  camera  data  and  sun  data, 
the  data  retrieval  components  for  selecting  appropriate  data  from  the  data  base, 
the  MERGE  program  for  merging  two  block  format  images,  and  the  CONV  program,  for 
converting  the  block  format  image  to  the  string  format  required  by  the  color  pass. 

All  of  these  components  are  under  the  control  of  an  executive  (PROGRAM 
TRAFFIC)  whose  function  is  to  read  and  respond  to  option  cards  and  thus  direct 
the  flow  of  a user  initiated  sequence  of  runs.  All  of  the  system  components 
mentioned  above  are  opted  by  the  user,  with  the  single  exception  of  the  data  re- 
trieval option,  for  which  the  executive  takes  the  responsibility.  The  user  also 
names  the  input  and  output  files  for  a run  and  thus  decides  whether  a current 
image  is  to  be  merged  with  a previously  produced  image  as  the  new  image  is  gen- 
erated or  whether  images  are  to  be  saved  separately  for  later  merges. 

The  only  card  input  required  by  the  system  are  the  option  cards  (two  for 
moat  options)  and  the  scene-setting  data;  camera  data,  sun  data,  translation  and 
scale  data  for  the  targets  and  nets,  raster  limits  for  each  geometry  and  shadow 
run,  and  color  assignment  cards  for  the  conversion  run.  For  a typical  whole 
scene  sequence  with  one  target  and  one  net,  there  will  be  four  geometry  runs, 
four  shadow  runs,  two  data  runs  and  one  conversion  run.  Apart  from  the  color 
assignment  cards  of  the  conversion  run,  fewer  than  fifty  cards  are  required. 
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Because  of  tHe  size  and  complexity  of  the  system,  it  has  been  necessary 
to  structure  it  in  overlays.  Since  the  overlay  method  is  machine-dependent, 
the  IBM  version  used  at  MAGI  in  the  production  phase  of  this  work  is  different 
from  that  required  by  to  the  government  for  the  CDC  6600  computer.  The  CDC 
version  is  the  subject  of  Section  5,  where  the  executive  and  new  routines  of 
the  system  are  described.  That  section  also  oontains  a user's  guide  and  input 
description. 

1.7  Demonstration  of  the  System 

The  completed  jystem  was  debugged  on  a data  base  consisting  of  a fifteen- 
tree  forest,  a tree  library  with  two  tree  types,  a terrain  based  on  nine  ele- 
vation points,  two  combinatorial  geometry  targets  and  a single  camouflage  net. 

A full  scale  problem,  based  on  a test  site  description  provided  by  the 
government,  was  then  executed  with  excellent  results.  The  details  of  that 
problem  are  described  in  Section  6 and  the  resulting  photograph  is  shown.  In 
that  section  we  also  give  execution  times  for  the  data-producing  codes  as  well 
as  the  image  generating  system,  and  make  recommendations  for  the  manner  in 
which  test-site  descriptions  should  be  provided  in  the  future  for  good  simu- 
lation. 


6 


2. 


THE  DATA  BASE 


2.1  Description  of  the  Data  Baae 

The  data  base  for  the  system  resides  on  a tape  or  disk  file  and  contains 
all  of  the  information  about  a single  test  site  as  well  as  all  of  the  data  for 
combinatorial  geometry  targets  and  camouflage  nets  that  may  be  used  in  connect- 
ion with  the  site. 

The  separate  libraries  or  data  units  which  congprise  the  data  base  are,  in 
the  order  of  occurrence  on  the  tape: 

a)  Forest  library 

b)  Tree  library 

c)  Terrain  library 

d)  Combinatorial  Geometry  target  library 

e)  Camouflage  Net  library 

The  detailed  format  of  the  data  base  is  given  in  Table  2.1.  Further  details 
of  the  forest  library  may  be  found  in  Section  3.5  of  this  report,  while  the 
terrain  library  is  fully  discussed  in  Section  4.2.3. 

The  forest  library  is  divorced  from  the  terrain  library  in  the  sense  that 
the  vertical  location  of  the  tree  base  is  at  sea  level.  The  actual  elevation  is 
extracted  from  the  terrain  library  at  the  time  of  retrieval  of  the  tree  from 
the  forest  library.  For  shadow  runs,  the  normal  to  the  terrain  is  also  computed 
at  this  time. 

Each  Combinatorial  Geometry  target  is  represented  by  card  images  of  the 
geometry  deck,  i.e.,  raw  data  prepared  by  the  user.  These  cards  are  fully  de- 
scribed in  references  (1)  and  (3).  A separator  record  consisting  of  the  Hollerith 
information  12X,FINX  folows  the  last  card  image  of  each  deck  except  the  last, 
which  is  followed  by  the  record  containing  the  information  12X,ST0P. 


Each  camouflage  rtt  requires  three  date  seta:  1)  a variable  number  of 
card  images,  2)  7 records  of  processed  binary  data  (the  output  of  CAM1IET  on 
TAPES , for  restart),  and  3)  two  final  card  images.  The  first  two  sets  are 
separated  by  a record  containing  the  Hollerith  information  12X,SWCH.  The  net 
library,  and  therefore,  the  data  base,  Is  closed  by  an  end  of  file. 

The  tree  library  describes  each  of  several  tree  types  by  means  of  two 
models:  a high  resolution  model  and  a statistical,  or  low  resolution  model. 

The  high  resolution  model  is  geometrically  identical  to  that  described  in  (1) 
and  (4)  but  much  of  the  multi-staging  has  been  removed  to  increase  the  efficiency 
of  ray-tracing.  The  low  resolution  model  is  based  on  the  high  resolution  model 
with  the  digitized  leaf  and  twig  structure  replaced  by  a statistical  leaf  cloud. 
These  models  were  prepared  by  the  government  and  are  described  in  Appendix  A 
of  this  report. 

The  identification  of  a model  ia  based  on  its  tree  type,  I(<_6),  and  its 
resolution.  Thus  the  identifier  for  the  high  resolution  modal  of  tree  type  X is 

LTYP (2*1-1)  - 10*1  + 1 } 

while  that  for  the  low  resolution  modal  of  the  same  tree  type  ia 

LTYP (2*1)  - 10*1  + 2 . 

This  identifier  is  the  first  record  of  a tree  model.  The  second  record  of 
the  model  carries  the  contents  of  the  labelled  coeaon  TREE  of  the  vegetation 
code  and  consists  of  construction  data.  A third  record  carries  some  statistical 
data  (meaningful  only  for  the  low  resolution  model),  namely  the  arrays  STTV(4), 
AVE  ( 3 ) and  DAVE  (2) , and  the  contents  of  the  cowaon  blocks  FG,  KD,KL0UD,  RAYTR 
and  COLR  of  that  model.  "here  are  21631  words  of  data  in  the  second  record 
and  1680  words  in  the  third  record. 
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There  are  currently  four  tree  types  represented  on  the  data  base:  a 

coniferous  tree  (LTYP(l)-ll,  LTYP (2)»12) , two  deciduous  trees  with  leaves,  the 
pin  oak  .(LTYP (3)-21,  LTYP(4)=22)  and  the  cottonwood  (LTYP(7)-41,  LTYP(8)=42), 
and,  finally,  a deciduous  tree  without  leaves  (LTYP(5)»31,  LTYP (6)“32) . Boom 
has  been  allowed  in  the  retrieval  routines  for  two  more  tree  types. 


TABLE  2.1  - FORMAT  OF  THE  DATA  BASE 


DATA  UNIT 

RECORD  NOS. 

CONTENTS 

comments 

1 

l.N 

FOREST  LIBRARY 

Output  of  FORGEN 

N+l 

78130 

Fla? 

2 

TREE  LIBRARY 

3 record*  for  each 

of  N models.  M^12 

1 

LTYP(l) 

Tree  ID 

2 

MA,FPD, .... 

3 

STTT, AVE. . . 

4 

LTYP (2) 

3*N-2 

• 

LTYP (N) 

Tree  ID 

3*N-1 

MA.FPD  .. 

3*N 

STTT, AVE. . 

3N+1 

15626 

Flag 

3 

1-9 

TERRAIN  LIBRARY 

Output  of  INSTAL 

10 

31252 

Flag 

4 

1 1 Nl 

CG  TARGET  LIBRARY 
Card  images  of 

NMAX  targets 

standard  input 

HOLLERITH  (20A4) 

Nl  + 1 

12X, FINIS 

Separator 

(Nl+2) , (NH-N2  + 2) 

Card  image* 

N1+N2+3 

• 

• 

12X, FINIS 
• 

• 

Separator 

• 

NMAX 

• 

E (N2+1) 
2-1 

12X,STOP 

Flag 

5 

NET  LIBRARY 

LMAX  nets 

1,L1 

LI  card  images 

for  restart  option  #5 

HOLLERITH  (20A4) 

Ll+1 

12X, SWCH 

Flag,  switch  to 
unformatted  data 

Ll+2 . Ll+8 

Processed  data 

7 binary  records 

from  CAMNET  restart  tape) 

L1+9,U+10 

2 card  images 

Final  cards  for  restart 

Repeat  above  for  each  net. 
End  data  baee. 
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2.2  Data  Retrieval  Methods  and  Routines 

The  data  retrieval  methods  ara  peculiar  to  tha  library  concerned.  Retrieval 
for  the  vegetation  modules  is  rather  complicated,  while  that  for  the  other  mod- 
ules is  essentially  a matter  of  tape  positioning. 

The  retrieval  routines  are  in  a primary  overlay  of  the  system,  OVERLAY (2), 
with  main  program  PROGRAM  DATRET.  This  overlay  is  called  by  the  executive  when- 
ever an  image  generating  or  shadowing  run  is  opted  by  tha  user.  The  geometry 
model  concerned  is  identified  by  a program  number,  JSAVE,  which  is  passed  to 
DATRET  through  common.  DATRET  calls  the  appropriate  subroutine: 

SUBROUTINE  TERRET 

Called  from  DATRET  for 

1)  terrain  image  or  shadow,  TERRET  searches  the  data 
base  for  the  integer  flag  15626  and  leaves  tha  tape 
positioned  for  reading  by  input  elements  of  the 
terrain  modules; 

2)  vegetation  data  processing,  TERRET  finds  the  terrain 
library  and  reads  into  memory  the  first  four  records 
(see  use  under  FORRET) . 

a)  Called  from:  DATRET 

b)  Subroutines  called:  None 

SUBROUTINE  TARGET (N,M) 

Called  from  DATRET  for  either  CG  targets  (M**l)  or  camouflage  nets  (M-2) , 
TARGET  searches  the  data  basa  for  the  appropriate  flag  and  positions  the  file  at 
tha  Nth  target  of  either  type.  The  data  are  later  read  by  the  input  routines 
of  the  appropriate  model.  The  ordinal  number  N is  supplied  by  the  user  as  de- 
scribed in  Section  5. 

s)  Called  from:  DATRET 

b)  Subroutines  called:  None 


SUBROUTINE  FORRET  * 

This  subroutine  is  callad  from  DATRET  to  retrieve  data  for  a particular 
model  of  a particular  tree  type.  If  the  model  is  not  in  the  tree  library  or  is 
not  represented  in  the  forest  library,  a flag  is  returned  to  the  executive 
inhibiting  the  geometry  run  and  calling  for  the  next  tree  model. 

The  action  of  FORRET  is  to  locate  the  model  on  the  data  base  and  copy  it 
to  a temporary  storage  device  (the  user  designated  logical  unit  IPD) . FORRET 
then  searches  the  forest  library  for  trees  of  the  correct  type  which  are  within 
the  scene  limits.  Each  such  tree  must  satisfy  a range  criterion  for  the  parti- 
cular model  (high  or  low  resolution) . To  calculate  the  range  and  to  complete 
the  tree  description,  the  vertical  position  of  its  bass  is  extracted  from  the 
terrain  library  by  a call  to  the  subroutine  ZANDW.  This  subroutine  will  also 
return  the  local  normal  to  the  terrain  if  a shadow  run  is  contemplated.  The 
twelve  or  (for  shadows)  fifteen  word  description  is  written  on  the  IPD  as  a 
single  record. 

The  range  criterion  for  high  or  low  resolution  model  is  specimen-dependent. 
For  a particular  tree  the  criterion  is  that  the  width  of  the  containing  box  shall 
subtend  at  the  observation  point  an  angle  greater  than  25  mils  to  rsquirs  the  use 
of  the  high  resolution  modal.  For  a tree  within  a dansaly  forested  region,  the 
value  50  mile  is  used.  Such  trees  ars  tagged  at  the  time  of  generation  as  de- 
scribed in  section  3.2. 

a)  Called  fromi  DATRKT 

b)  Subroutines  calledt  ZANDW, INROUT 
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SUBROUTINE  ZANDW  . 

This  is  the  routine  called  by  FORRET  to  calculate  the  elevation,  Z,  of  the 
terrain  at  the  position  (x,y)  of  a tree  base.  For  shadow  runs,  ZANDW  will  also 
find  the  local  normal  to  the  terrain.  For  these  calculations,  ZANDW  calls 
routines  of  the  terrain  ray  tracing  module. 

a)  Called  fromi  FORRET 

b)  Subroutines  called:  INXPOS,  NORM 

SUBROUTINE  INROUT (JFL) 

This  subroutine  is  called  by  FORRET  to  determine  for  a particular  group  of 
trees  contained  within  known  rectangular  limits  (see  Section  3),  whether  any 
part  of  the  rectangle  is  within  the  scene  limits.  The  flag  JFL  is  returned  as 
non-zero  if  the  answer  is  positive. 

a)  Called  from:  FORRET 

b)  Subroutines  called:  None 

2 . 3 Forest  Data  Organising  Programs 

Because  of  the  very  large  memory  requirements  for  a tree  model  (over  23000 
words) , it  has  been  neces3ary  to  design  the  vegetation  runs  around  a single  model 
and  to  limit  the  number  of  trees  in  memory  at  one  time  to  two  hundred.  A complete 
ground  cover  image  is  obtained  by  a loop  on  model  number,  with  a continuous  merge 
of  the  image  in  process  with  that  of  the  previously  processed  models.  This 
"merge-on-the-fly"  process  is  described  in  Section  5.1.  The  loop  extends  from 
JFIRST  (input  or  default,  1,  to  JLAST,  input  or  default,  12). 

The  limit  on  the  number  of  trees  concurrently  in  memory  is  observed  by 
organizing  the  trees  of  the  scene  by  blick  number  of  the  image  scan  (line  by 
line,  left  to  right  within  a line).  This  function  is  performed  by  the  PROGRAM 
PREPAR  for  image  generation  and  by  the  PROGRAM  PREPR  for  shadowing  by  the  forest. 


The  organized  treefi  are  written  on  TAPE3  in  records  of  200  trees,  each  record 
headed  by  an  identifying  block  number.  Clearly  trees  may  be  repeated  for 
blocks  qf  different  number.  A second  output  tape  of  PREPAR,  TAFE10,  carries 
for  each  covered  block  an  array  of  dimension  2 by  20,  one  column  of  which  orders, 
on  distance  from  the  camera,  the  tree  numbers  to  be  investigated  for  that  block. 
The  other  column  gives  the  distance. 

These  programs  operate  on  the  TAPE  IPD,  the  output  of  FORRET,  which  con- 
tains those  specimens  of  a single  model  which  are  within  the  scene  limits  estab- 
lished by  the  subroutine  COVER  (scene  setting  overlay). 
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3. 


THE  FOREST  GENERATOR:  PROGRAM  FORGEN 


3.1  Introduction 

The  purpose  of  this  stand-alone  program  is  to  create  the  data  file  which 
contains  the  type,  size,  location  and  orientation  for  every  tree  in  the  area 
of  interest.  FORGEN  performs  this  function  in  response  to  relatively  simple 
input  data  by  selecting  values  from  internal  distribution  functions  with  user- 
defined  means  and,  in  some  cases,  user-defined  variances.  The  output  data  are 
organized  by  region  (defined  below)  and,  within  each  region,  by  tree  type.  The 
data  for  an  individual  tree  are  given  as  the  parameters  of  the  containing  box, 
i.e.,  the  vector  position  of  one  base  corner  and  three  vectors  specifying  the 
edges  of  the  box,  and  therefore  its  size. 

In  the  following  sub-sections  of  this  section,  we  will  describe  the  phil- 
osophy and  design  of  the  program,  the  subroutines  comprising  the  program,  the 
input  preparation  and  logical  unit  assignments  and  the  format  of  the  output  file. 

3.2  Program  Design 

The  design  of  the  program  was  guided  by  observations  on  the  kinds  of  tree 
distributions  that  might  be  of  interest  and  by  practical  considerations  such  as 
the  possibility  of  replacing  individual  parts  of  the  description  or  adding  to  the 
description  without  rerunning  the  entire  problem. 

Three  kinds  of  tree  configurations  were  identified  as  useful  to  efficient 
generation  of  realistic  ground  cover.  A basic  concept  in  the  description  of 
two  of  these  configurations  is  the  "cluster"  and  is  based  on  the  fact  that,  in 
low  density  distributions  of  trees,  the  unit  grouping  may  contain  more  than  one 
tree.  Thus,  for  our  purposes,  trees  may  occur  in  clusters  whose  size  is  de- 
termined from  an  input  mean  and  a resident  distribution  function.  The  three 


configurations  are: 


1. 


A uniform  distribution  of  clusters  of  each  of  several  types  of 


trees  within  a cl6sed  boundary.  The  user  defines  the  boundary  and  the  types 
of  trees  and,  for  each  tree  type,  the  mean  cluster  size  (up  tc  10),  the 
fractional  density,  the  mean  height  and  the  half-width  of  the  height  range. 

2.  A uniform  distribution  of  clusters  of  each  of  several  tree  types 
along  a "thick"  line.  The  user  defines  the  line,  its  thickness  and  tree  para- 
meters as  in  (1)  above. 

3.  A distribution  of  individual  trees  whose  types»  sizes,  locations  and 
some  aspects  of  orientation  are  defined  by  the  user. 

The  first  of  these  clearly  lends  itself  to  the  description  of  forested 
areas  or  fields  with  scattered  trees  or  tree  clusters.  The  second  might  repre- 
sent the  fringe  arrangement  that  frequently  occurs  near  a road  or  on  the  bound- 
aries between  fields  in  farm  country.  The  third  type  of  representation  is  use- 
ful in  the  immediate  neighbornood  of  a target,  since  a very  special  local  con- 
figuration may  be  selected  for  camouflage  purposes. 

In  what  follows,  we  will  distinguish  these  configurations  by  their  geometric 
peculiarities  and  refer  to  them  as  "boundaries",  "lines",  and  "point  sets”,  re- 
spectively. The  entire  organization  of  the  program  - input,  program  action  and 
output,  as  well  as  the  replacement  and  addition  of  data  - is  based  on  this  three- 
part  description.  The  order  of  input,  treatment  and  output  is  always:  boundaries, 
lines,  point  sets. 

Because  the  line  segments  comprising  the  boundary  and  the  line  may  be  deter- 
mined by  a road  or  a stream  bed,  the  user  may  give  a (positive)  displacement  for 
each  segment  such  that  the  segment  used  by  the  codes  is  clear  of  the  road  or 
stream  bed.  Thus  the  displacement  will  normally  be  somewhat  more  than  half  the 
width  of  the  bed.  The  direction  of  displacement  for  a boundary,  whose  defining 
points  must  be  given  counter-clockwise,  is  always  inward.  The  line  must  be  de- 


scribed such  that  the  directed  segment  from  point  N to  point  N+l  is  displaced 
to  its  left.  These  ideas  ars  illustrated  in  Figure  3.1. 


A restriction  on  the  closed  boundary  should  be  mentioned  here,  namely 
that  a line  X = .constant  should  cut  the  boundary  in  two  ard  only  two  points. 

If  a bounded  region  does  not  satisfy  this  restriction,  it  should  be  sectioned 
into  two  or  more  bounded  regions. 

The  distribution  of  trees  within  a boundary  proceeds  by  determining  the 
smallest  cont a Lning  rectangle  with  sides  parallel  to  the  x-y  axes  and  covering 
this  area  with  a square  grid.  The  grid  size  ia  the  product  of  the  average 
width  of  a tree  box  in  the  area  by  an  input  factor  which,  in  effect,  determines 
the  mean  spacing  of  trees  within  a cluster.  If  a tree  is  assigned  to  a ceil 
of  the  grid  the  base  location  will  be  chosen  somewhere  within  the  cell. 

Once  the  grid  has  been  established,  an  array  is  set  up  which  establishes, 
for  each  value  of  the  x index,  I,  the  limits  of  the  y index,  J,  that  correspond 
to  the  boundary.  Simultaneously,  a pointer  array  is  set  up,  such  that  for 
each  value  of  I,  the  emulative  number  of  grid  points  through  1-1  is  recorded. 
These  pointers  refer  to  the  "IQ"  array,  which  will  contain  the  identifiers  for 
trees  assigned  to  grid  points. 

The  assignment  of  trees  to  calls  of  the  grid  requires  that  cluster  den- 
sities be  established  for  each  tree  type.  The  cluster  density  is  sampled  for 
a grid  point  which  will  act  as  an  "anchor"  pcint  for  location  of  a cluster. 

The  cluster  size  is  then  determined  from  the  cluster  mean  and  a highly  peaked 
distribution  function  which  covers  a range  from  unity  to  twice  the  mean.  A 
square  array  of  grid  points  is  than  investigated  for  unoccupied  sites  and  the 
members  of  the  cluster  are  assigned  to  these  with  probability  that  decreases 
with  distance  from  the  center.  This  has  the  effect  of  softening  the  corners. 

It  is  possible  that  some  members  of  the  cluster  will  not  be  accommodated  in 
the  array,  but  this  loss  is  compensated  at  the  end  of  the  assignment  process 
by  distributing  individual  trees.  Tha  dimension  of  the  square  array  is  de- 
termined from  the  expression 

LENGTH-SQRT (2* (KL+2) ) 


1 1 


FIGURE  3.1  - Boundaries  and  "Thick"  Lines 
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(a)  A closed  boundary:  points  1 through  5 lie  along  a 
road  bed.  Non-zero  values  of  W(l)  through  W(4) 
give  the  appropriate  displacements. 


Input  lino 


(b)  A line:  points  1 through  5 lie  along  the  same  road  bed 
aa  in  (a) . By  describing  these  segments  in  reverse,  tha 
displacement  is  forced  in  the  opposite  direction. 


where  KL  is  the  clyster  size.  If  the  mean  cluster  size  is  unity  this  method 
is  by-passed  and  an  individual  tree  is  located  at  the  selected  grid  site. 

A similar  treatment  is  applied  to  each  segment  of  a "thick"  line,  with 
the  exception  that  the  members  of  a cluster  occupy  successive  grid  sites  along 
the  line,  each  site  selected  from  a stack  which  is  perpendicular  to  the  line 
and  gives  the  lines  its  "thickness".  Thus  the  thickness  is  used  to  "scatter" 
the  trees  perpendicular  to  the  line. 

The  fractional  density  given  by  the  user  for  each  tree  type  can  now  be  de- 
fined* for  the  bounded  region  it  ir  the  fraction  of  grid  sites  that  will  be 
occupied  by  this  type?  for  the  line  it  is  the  fractional  number  of  stacks  which 
is  occupied  by  treeB  of  this  type,  and  one  o a stack.  If  a line  is  composed  of 
trees  of  more  than  one  type  it  is  possible  that  a stack  will  have  more  than  one 
occupant. 

The  size  of  a tree  is  chosen  from  a flat  distribution  abou‘  the  given  mean 
and  of  half-width  given  by  the  user.  The  polar  angle  of  the  trunk  is  chosen 
uniformly  within  a cone  whose  ha  If -angle  is  0.1  radian.  Both  U.e  azimuth  cf  the 
trunk  and  the  axial  rotation  about  the  trunk  are  chosen  from  uniform  distributions 
between  0.0  and  2ir.  The  location  of  the  center  of  the  base  of  the  containing  nox 
is  selected  uniformly  within  a circle  about  the  grid  site  with  radius  equal  to 
half  the  grid  spacing. 

Each  tree  associated  with  a point  set  is  fully  described  by  the  user  except 
for  the  axial  rotation  about  the  trunk.  This  latter  parameter  is  selected  as 
for  trees  in  the  other  configurations.  However,  it  was  felt  that  a user  might 
wish  to  describe  a fallen  or  leaning  tree  in  soma  specific  azimuthal  direction 
and  specification  of  the  polar  and  azimuthal  angles  is  necessary  for  this.  If 
these  entries  are  left  blank  the  trees  will  be  upright  but  the  axial  rotations 
will  give  aufficient  variation  in  appearance. 
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A rough  printer  plot  capability  is  provided  within  the  program  and  is 
applied  to  each  bounded  region,  each  segment  of  a line  and  each  point  set.  For 
the  bounded  region  and  for  the  line  segment  (always  represented  horizontally) 
the  technique  is  to  identify  grid  points  with  printer  positions  and  for  each 
occupied  grid  site  to  print  the  tree  identifier  which  is  an  integer  between  1 
and  6.  The  same  procedure  is  used  for  representing  the  point  set  except  that 
the  grid  size  is  supplied  by  the  user. 

A better  plotting  capability  is  provided  outside  the  program  by  the  PROGRAM 
PLOTFR.  The  card  input  u the  program  ia  a single  card  giving  the  rectangular 
limits  of  the  plot  and  a grid  size.  PLOTFR  operates  on  the  output  tape  of  FORGEN 
to  provide  an  integrated  view  of  the  specified  region. 

The  identification  of  a tree  on  the  output  tape  of  FORGEN  do.  j not,  of 
course,  distinguish  between  the  high  resolution  model  and  the  low  resolution 
model  for  that  tree.  That  determination  must  be  made  for  each  siene,  i.e.,  for 
each  ratio  of  tree-size  to  distance  from  the  camera.  However,  t •»  criterion  for 
switching  from  one  model  to  the  other  may  properly  be  influenced  by  the  tre“ 
situation.  For  example,  a tree  in  the  interior  of  a dense  forest  need  not  be 
represented  as  carefully  as  a single  specimen.  With  this  in  mind,  FORGEM  tags 
such  trees  by  attaching  a minus  sign  to  an  intrinsically  positive  quantity, 
namely,  the  vertical  component  of  the  height  of  the  containing  box.  The  basis 
of  this  assignment  is  that  the  tree  must  be  a non-boundary  tree  of  a bounded 
region  for  which  thn  fractional  area  occupied  by  trees  is  greater  ~nan  0.2. 

3.3  Description  of  Routines 

In  this  section,  we  describe  all  of  the  routines  of  FORGEN.  A linkage  chart 
is  given  in  Figure  3.2. 
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PROGRAM  FORGEN 


FORGEN  governs  the  generation  of  the  forest  library  and  performs  many  of 
the  required  calculations.  In  particular,  it  distinguishes  between  the  gener- 
ation of  a new  library  and  the  updating  of  an  old  library,  and  reflects  the  organ- 
ization of  the  calculation  according  to  boundaries,  lines  and  points. 

a)  Main  program 

b)  Subroutines  called:  DISPL,  INPROC,  GRID,  ASSIGN, 

TCAL,  BOXIT,  POSIT,  PLOTl 

SUBROUTINE  ASSIGN 

ASSIGN  is  called  from  the  main  program  for  boundaries  and  thick  line  seg- 
ments after  the  grid  has  been  established.  It*s  function  is  to  select  a cluster 
size  and  a cluster  position  and  to  assign  the  members  of  the  cluster  to  grid 
sites.  ASSIGN  does  this  for  the  number  of  clusters  required  to  provide  the  cal- 
culated cluster  density.  Because  the  number  of  necessary  sites  are  not  always 
available  in  the  cluster,  ASSIGN  finds  individual  sites  to  complete  the  assign- 
ment. Since  the  assignment  process  is  different  for  boundaries  and  line  seg- 
ments, a parameter  NOPT  {>0  for  boundaries,  *1  for  line  segments)  is  passed  to 
ASSIGN  from  the  main  program. 

a)  Called  from:  FORGEN 

b)  Subroutines  called:  CLSTR 

SUBROUTINE  BOXIT 

This  subroutine  calculates  the  box  parameters,  i.a.,  the  vertex  vector  and 
three  edge  vectors,  for  an  individual  tree,  assuming  an  elevation  of  z»0.0  for 
the  base.  It  selects  an  axial  rotation  angle  and  accepts  from  the  calling 
routine  a horizontal  base  position,  a polar  angle,  an  azimuth  and  a magnification 
of  the  containing  box  for  that  tree  type. 

a)  Call sd  from:  FORGEN,  TCAL 

b)  Subroutines  called:  None 
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SUBROUTINE  CLSTR 

CLSTR  selects  a cluster  size  and  returns  thin  value  to  ASSIGN.  The  dis- 
tribution function  is  symmetric  about  the  mean  with  range  from  unity  to  one 
less  than  twice  the  mean.  The  probabilities  decrease  by  a factor  of  one  half 
with  each  additional  unit  departure  from  the  mean. 

a)  Called  fromt  ASSIGN 

b)  Subroutines  called:  None 

SUBROUTINE  GRID 

It  is  the  function  of  GRID  to  set  up  {for  bounded  regions  only)  the  array 
which  defines  the  bounded  region  in  terms  of  grid  indices.  LIST(1,I)  is  the 
lowest  index  (Jl)  for  horizontal  grid  lines  passing  through  the  vertical  grid 
line  of  index  I,  while  LIST(2,I)  is  the  highest  such  index  (J2) . LIST(3,I)  is 
the  cumulative  sum  of  (J2-J1+1)  through  1-1  and  thus  constitutes  a pointer  to 
an  array  carrying  the  grid  site  assignments.  The  process  is  so  trivial  for 
thick  line  segments  and  for  point  sets  that  it  is  performed  in  FORGEN. 

a)  Called  from:  FORGEN 

b)  Subroutines  called:  Nona 

SUBROUTINE  INPROC 

This  subroutine  ie  celled  from  FORGEN  to  read  the  data  for  all  the  trees 
of  a bounded  region  or  a line,  i.e.,  number  of  types,  type  identifier,  mean  cluster 
size,  fractional  density,  mean  height  of  containing  box,  relative  range  of  heights 
for  this  region,  and  the  siurink  factor.  It  than  calculates  the  cluster  prob- 
abilities mean  magnifications  (ratio  of  input  hsight  to  resident  box  height  for 
that  tree  type)  and  the  grid  size. 

a)  Called  from:  FORGEN 

b)  Subroutines  called:  None 
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SUBROUTINE  FLOT1 

PL0T1  is  used  to  display  on  the  printer  the  grid  site  assignments  for 
each  region.  The  position  in  a printer  line  corresponds  exactly  to  the  x-index 
I,  while  the  line  itself  corresponds  to  the  y-index  J.  The  grid  spacing  is 
calculated  for  the  boundary  or  line  segment  but  is  an  input  for  the  point  set. 
PL0T1  operates  from  the  arrays  LIST  and  IQ.  LIST  has  already  been  described. 
IQ(K)  contains  a zero  for  an  unoccupied  site  and  an  identifying  integer  from 
1 to  6 for  an  occupied  site.  LIST{3,I)  is  a value  of  K such  that  K+J  is  the 
location  in  IQ  corresponding  to  the  (X,J)  grid  site.  Segments  of  lines  are 
plotted  horizontally,  but  actual  end  points  are  printed  out. 

a)  Called  froms  FORGEN 

b)  Subroutines  called:  None 
SUBROUTINE  POSIT 

This  routine  is  called  by  FORGEN  whan  the  code  is  being  operated  in  the  up- 
dating mode,  i.e.,  the  mode  in  which  an  existing  library  is  to  have  certain  ele- 
ments replaced  and  other  elements  added.  POSIT  must  copy  acceptable  elements 
from  the  old  tape  to  tne  new  tape  arid  position  the  old  tape  past  unacceptable 
elements. 

a)  Called  from:  FORGEN 

b)  Subroutines  called:  None. 

SUBROUTINE  TCAL 

TCAL  governs  the  calculation  of  tree  box  parameters  for  boundaries  and 
lines.  It  selects  the  final  location  of  ths  bass  of  ths  tres  within  the  grid 
cell,  for  all  trees  of  ray  regiona,  a magnification  within  the  input  limits  about 
the  mean  magnification,  a polar  angle  and  an  azimuth.  BOXIT  is  then  called  to 
J complete  the  calculation  for  one  tree.  TCAL  also  writes  out  on  tape  (TAPE2)  the 

12-word  tree  descriptions  in  records  of  200  trsss.  Ths  calculation  is  ordered  on 
/ tree  types. 

' a)  Called  from:  FORGEN 

j 

b)  Subroutines  called:  BOXIT 


i 
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3.4  Input  and  Logical  Unit  As^lg7unent3 

All  input  data  are  read  from  cards  on  one  of  two  formats: 
# 

Integer  data:  FOIiHAT(dllO) 

Floating  Point  data:  format  (6E12. 4) 


CARD  NO.' 


CONTENTS 


COMMENT 


1 


ISUBST 


NB 
NL  > 

NP 

J 


Mode  option  parameter  equal  to  0 or  1. 

0:  create  new  library  from  scratch; 

1:  update  existing  library. 

Number  of  (boundaries,  lines,  point  sets) 
to  be  processed  for  ISUBST-0;  number  to 
be  replaced  for  ISUBST*!. 


Cards  2,3  and  4 are 


included  only  for  I6UBST**!. 


2 

3 

4 


IDB(I) ,I-1,NB+1 
IDL(I) ,I-1,NL+1 
IDP  (I)  , I*1,NP+1 


I^(NB,NL,NP) , ordinal  zv  nber 
of  (boundary,  line,  point  set)  to  be 
replaced.  I» (NB+1,NL+1,NP+1)  number  of 
such  regions  to  be  added  (may  be  zero) . 


Card  sequence  S through  11  is  needed  for  each  bounded  region  end  line  with  the 


exception  of  Card  7 which  ia  needed  only  for  lines.  Card  sequence  10,  11  is  re- 


peated for  each  tree  type  in  each  description. 


5 

NXY 

% 

Number  of  points  used  to  describe  this 
boundary  or  line. 

6 

XPl 

YP) 

X and  Y coordinates  of  a point; 

w 

displacement  of  line  segment  between  that 
point  and  the  next.  TVo  points  per  card. 

7 

WN 

for  lines  only,  "thickness"  or  allowed 
straggle  of  bases. 

8 

NTYP 

Number  of  tree  types  for  this  region. 

9 

SHRINK 

Factor  to  be  applied  to  tree  box  width  to 
determine  swan  spacing  of  trees  in  a cluster 

10 

ITYP 

'XMP 

Tree  identifier  (11,21,31,41,51,61); 
moan  cluster  size  for  this  type;  one  type 
por  card. 

11 

FRCT 

Fraction  of  grid  points  covered  by  this 
tree  type; 

SUN 

SIC 

mean  height  of  this  tree  type;  t actional 
height  variation  for  this  tree  type.  One 
type  par  card. 

as 


Card  sequence  12  through  15  is  needed  for  each  point  set.  Card  15  is  repeated 
for  each  tree  in  a point  set  with  order  of  tree  type  given  by  card  13. 


CARD  NO 

CONTENTS 

COMMENT 

12 

NXY 

Number  of  trees  in  this  set,  number 

NTYP 

of  tree  types. 

13 

ITYP 

Tree  type  identifier; 

NUMR 

number  of  trees  of  this  type.  Three 
sets  per  card. 

14 

xminI 

XMAX\ 

x-limits; 

yminL 

YMAX) 

y-limits; 

DS 

grid  spacing 

15 

*1 

x,y  coordinates  of  tree 

base; 

SMN 

height  of  this  tree; 

THET 

polar  angla  of  cfam  (rad? ; 

PHI 

azimuth  of  stem  (rad). 

A final  note  on  input  cards  5 and  6 for  boundaries:  the  closing  is  produced 
by  the  code  on  the  assumption  that  the  last  point  is  not  a repetition  of  the  first 
point.  For  both  boundaries  and  lines,  the  quantities  W(I)  should  be  made  some- 
what greater  than  half  the  width  of  the  road  or  stream  bed  since  it  is  the  base 
of  the  tree  which  is  confined  to  be  within  the  boundary  or  within  the  "thickness" 
of  tho  line.  Finally,  the  points  defining  a boundary  must  be  sequential  and  in 
counter-clockwise  order,  while  those  defining  a line  must  be  described  sequentially 
in  a direction  such  that  the  stream  or  road  bed  is  on  the  right.  The  w values 
are  all  positive. 

Logical  units: 

TAPE1  : existing  tree  library  for  update  only  (ISUBST-l) 

TAPES  i output  tree  library. 
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3.5  The  Output  Tape 

The  output  tape  (logical  unit  2)  is  organist  tf  region  typer  bounded 
regions,  lines  end  point  sets.  Within  each  type  the  order  is  that  determined 
on  input. 

RECORD  NO.  CONTENTS  COWttjTg 

1 NBT , NLT , NPT  Numbers  of  Boundaries,  lines  and  point  sets 

Each  region  is  represented  by  the  following  sequence  of  records. 

2 IDUM,NTYP,  IDUM2  f/TYP  - number  of  tree  types 

( ITYT  ( I ) , NUMR  ( X ) ,1-1, NTYP ) Type  identifier,  number  for  all  types 

XMIN,XMAX,YMIN, YMAX  Rectangular  limits  of  region 

Each  tree  type  in  a region  ia  represented  by  the  following  sequence  of  records  t 

3 (TT(I,J)  ,1-1,12)  ,J-1, 200)  Up  to  200  tree  boxes  for  first  tree 

„ type  in  each  record. 

e 

m 

L1-2*{NUMR(1)/20(HU 

Each  tree  type  starts  a new  record. 
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4 . THE  TERRAIN  PROGRAM 

4.1  Introduction 

The  terrain  model  that  was  developed  under  a previous  contract^-  was  a 

2 

modification  of  the  camouflage  net  model  that  embodied  softened  contour  lines 
and  a "road"  capability.  It  was,  however,  wedded  to  inflexible  color  pattern 
boundaries,  characteristic  of  the  nets,  and  to  the  input  mode  of  the  net,  which 
is  quite  inappropriate  to  a terrain  model.  The  purpose  of  the  effort  reported 
here  was  to  remove  these  restrictions  from  the  model,  to  separate  the  establish- 
ment of  the  terrain  library  from  the  ray  tracing  process,  and  to  provide  for  easy 
retrieval  of  the  data  for  ray  tracing.  The  first  two  steps  required  the  develop- 
ment of  the  stand-alone  program,  PROGRAM  INSTAL,  and  the  replacement  of  a single 
routine  in  the  ray  tracing  module.  The  last  step  is  a matter  of  positioning  the 
data  base  tape  and  reading  the  processed  data  into  memory. 

A further  deficiency  of  the  old  model  was  an  extremely  slow  trianglulation 
method  which,  in  the  course  of  establishing  the  reel  terrain  of  nearby  one 
hundred  points,  became  intolerable.  Modification  of  the  subroutine  SEGNET  was 
necessary  to  overcame  this  problem. 

4.2  PROGRi-M  INSTAL 

4.2.1  Design 

This  is  a stand-alone  program  for  reading  and  processing  elevation  data, 

* 

road  data,  color  boundary  data  and  region  number  assignments.  The  output  is  a 
file  which  is  to  constitute  the  third  library  on  the  data  base  tape.  The  file 
consists  essentially  of  the  common  blocks  required  by  tha  ray  tracing  routines. 

In  the  previous  model,  the  terrain  installation  procsss  took  place  in  two 
steps*  1)  establishment  of  the  flat  "net"  description  in  2-space;  2)  construct- 
ion cf  the  installed  "net"  in  3-space  and  calculation  of  parameters  for  mapping 
a (strike)  point  on  the  installed  net  onto  a point  of  the  flat  net  for  (struck) 
region  identification.  PROGRAM  INSTAL  replaces  the  first  step  entirely  end 
modifies  same  portions  of  the  second  step  so  that  simple  elevation  data  ere  used 
in  the  installation  process. 
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The  "flat  net"  is  essentially  replaced  by  supplying  a projected  region 
boundary,  i.e.,  a sequence  of  points  (x,y)  in  a horizontal  plane  defining  a 
closed  boundary  on  the  elevated  terrain.  Any  terrain  not  enclosed  by  such  a 
boundary  is  lumped  under  a single  region  assignment.  In  the  ray  tracing  pro- 
cedure, the  "mapping"  problem  is  replaced  by  the  problem  of  determining,  for 
given  (x,y),  the  corresponding  boundary  number. 

The  treatment  of  roads  has  not  been  changed  except  that  the  road  input 
reading  and  processing  is  now  part  of  PROGRAM  IN3TAL. 

4.2.2  Description  of  Routines 

Table  4.1  lists  the  routines  of  PROGRAM  INSTAL  and  identifies  them  as  new, 
modified  or  unchanged.  Only  the  first  two  categories  will  be  described  here. 

PROGRAM  INSTAL 

This  is  the  main  program  whose  function  is  to  call  the  specific  input  read- 
ing and  processing  routines  and  to  write  the  processed  data  (essentially  the 
c canon  blocks)  out  onto  TAPES. 

a)  Main  program 

b)  Subroutines  culledi  IHPROS,  S EG NET,  TRINET,  FLIP,  NETBOD, 

ROADIN,  INBOUN 

SUBROUTINE  INPROS 

This  routine  reads  elevation  data  (x,y,z)  for  ISUP  interior  points,  and 
locations  (x,y)  of  NGM  sea  level  (mO.O)  points  on  an  enclosing  boundary.  The 
output  in  the  HPT  array  containing,  for  up  to  100  points,  the  five  word  greup 
x,y,s,u,v  where  u-x  and  v-y . This  redundancy,  produced  by  replacing  the  flat 
net  description,  does  no  harm  and  allows  us  to  retain  the  other  processing  and 
ray  tracing  routine*  intact. 
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Because  elevation  values  (z)  are  usually  provided  in  feet,  while  hori- 
zontal locations  and  distances  are  given  in  meters,  the  code  assumes  this 
mixed  unit  input  and  divides  the  elevations  by  3.048  ft/meter, 

a)  Called  from:  INSTAL 

b)  Subroutines  called:  None 

SUBROUTINE  SECNET 

This  routine  receives  the  HPT  array  from  INPRQS  and  performs  the  triangu 
lation  of  the  terrain.  The  method  is  to  form  all  possible  segments  between 
pairs  of  interior  and  to  order  the  segments  according  to  increasing 

length  for  each  nornt.  The  shortest  segment  is  accepted  as  a triangle  side 
if  it  does  not  cross  any  previously  accepted  segment  and  is  not  parallel  to 
a previously  determined  side  of  the  current  triangle.  The  necessary  improve- 
ment introduced  into  this  routine  is  a sort  on  segment  length  prior  to  the 
point  by  point  search  for  triangle  sides. 

a)  Called  from:  INSTAL 


b) 


Subroutines  called:  None 


TABLE  4.1  - The  RouUms  of  PROGRAM  INSTAL 


NAME 

COMMENTS 

INSTAL 

Main  prog  ram  % new 

INFROS 

Modified 

SEGNET 

Modified 

TRINET 

Old 

FLIP 

Old 

NETBOD 

Old 

ROAD IN 

Old 

INBOUN 

Hew 

CROSS 

Old 

♦:orm 

Old 

INIPOS 

Old 

XYZ2UV 


Old 
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SUBROUTINE  INBOUN 


This  routine  treads  and  process**  region  boundary  data  and  makes  region 


number  assignments  according  to  input.  A sequence  of  NXY  points,  (x,y)  counter- 
clockwise about  a closed  boundary  (bad  point  not  repeated)  is  read  for  each  of 
NB  boundaries.  A pair  of  direction  cosines  is  computed  and  stored  for  each  line 
segment.  The  rectangular  limite  of  the  boundary  are  also  stored  to  speed  up  the 
search  process  which  takes  place  in  ray-tracing.  A list  of  HB+1  color  region 
numbers  is  read  and  stored,  the  last  corresponding  to  all  uncovered  terrain. 

a)  Called  from*  INSTAL 

b)  Subroutines  called:  None 

4.2.3  Input,  Output  and  logical  Unit  Assignments 


Card  (1)  LStIP,  NUM 


F0RMAT(2I5) 


LSUP 

NUM 

Card  (2)  x,y, z 


Card  (3)  x,y 


Number  of  interior  points  for  elevation  data 

Number  of  boundary  points  with  z<=0.0. 

Read  by  INPROS. 

FORMAT (3F10.0) 

Elevation  data  for  bSUP  points,  1 point  per  card. 
Read  by  INPROS.  Scale  difference  between  (x,y) 
and  z assumed:  (x,y)  in  meters,  z in  feet. 

FORMAT (2F10.0) 

Locations  of  NUM  boundary  points,  1 point  per  card. 
Read  by  INPROS. 


Cards  (4)  through  (6)  Read  by  HOADIN  - see  (1)  p.  43. 

Card  (7)  NB  FORMAT (15) 

Number  of  region  boundaries  to  be  read. 

Read  by  INBOUN. 

Card  sequence  (8)  and  (9),  read  by  INBOUN,  repeated  for  NB  boundaries. 


Card  (8)  NXY  F0RMAT(I5) 

Number  of  points  for  this  boundary 


Card  (9)  XP,YP  FORMA7(2F10„0) 

Coordinates  of  NXY  points,  1 point  per  card. 

Card  (10)  ICOLOR(I)  ,I-1,NB+1  FORMATU4I5) 

Ordered  color  numbers  for  NB  bounded  regions  and 
"rest  of  world".  14  per  card.  Read  by  INBOUN. 
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Logical  unit  assignment:  TAPE8  «=  output  tape. 

The  output  tape  has  nine  records,  each  of  which  contains  the  data  for  a 
labeled  common  block  or  data  for  unlabeled  comnon  as  used  in  the  ray  tracing 
elements  of  the  model.  The  records  are  as  follows: 


RECORD  NO. 

CONTENTS 

COMMON 

1 

LSUP,IPEG,LIHT,HPT(500) 

LIMAGE 

. 

2 

LIJ (600) 

HSEGNE 

3 

NT, LE TRIP (600) 

HTRINE 

4 

/V, HBD (4000) 

HTKXNE, 

unlabeled 

comnon 

5 

jit,c.tit,  lrn 

OOMRN 

6 

NUMROD , NHIT , NCROAD 

ROADC1 

7 

IPROAD (25) , IPNORM (25) ,NSEG2(25) 
WIDTH (25) , WIDTHS (25) ,NCR(25), 
LRUT ( 25) . FRAC (25) ,CPAC(25) 

ROADC2 

8 

xp (looo) ,yp (looo) ,wx (looo) , 

WY(iCOO) ,NP(100) ,NB 

POINT 

9 

1J 

4.3  Terrain 

XLIST(400) ,ICOLOR(100) 
EOF 

Ray-Tracing 

ORGAN 

The  ray-tracing  module  of  the  new  terrain  model  differs  front  that  of  the 
old  model,  described  in  a and  b,  only  in  the  replacement  of  the  routine  FIND 
by  a pair  of  routines,  TFIND  and  LOCAL  (called  by  TFIND) . 
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SUBROUTINE  TFIND  . 

This  routine  is  called  from  RAYTRK  after  a ray  has  strucK  the  terrain  at 
a point  (x, y,z).  Using  the  coordinate  x and  y,  FIND  examines  the  rectangular 
limits  of  each  of  the  bounded  regions  to  determine  whether  that  region  may 
contain  the  point.  If  a region  is  a candidate,  FIND  calls  LOCAL  to  perform  the 
search  with  the  real  boundary.  If  a bounded  region  IR  contains  the  point,  the 
assigned  region  number,  ICOLOR(IR)  is  extracted  and  returned  to  RAYTRX.  'Other- 
wise it  is  assumed  that  the  strike  point  is  in  the  "rest  of  the  world"  and 
ICOLOR(NB+l)  is  returned. 

a)  Called  from:  RAYTRK 

b)  Subroutines  called:  LOCAL 

SUBROUTINE  LOCAL (IR,IFL) 

This  routine  uses  the  ordered  sequence  of  points,  (X,P(I) ,YP(I) ,I=1,NXY) , 
that  define  the  boundary  of  the  region  IR  and  the  associated  direction  cosines, 
(WX (I) ,WY  (I) ) , to  determine  whether  the  point  x,y  is  within  the  region.  The 
number  of  boundary  points  NXY  is  extracted  from  the  array  of  pointers  NP.  If 
the  point  is  within  the  region  the  signal  IFL  is  set  to  unity?  otherwise  a zero 
is  returned  to  FIND. 

a)  Called  from:  FIND 

b)  Subroutines  called:  None. 

4 . 4 Data  Retrieval  for  Ray  Tracing 

* 

The  data  retrieval  for  the  terrain  model  is  under  the  control  of  the 
executive  of  the  overlay  system  and  is  discussed  in  Section  5.  Since  the  pro- 
cessed terrain  data  file"  is  the  third  date  library  on  the  1DB,  retrieval  con- 
sists in  positioning  the  IDB  (SUBROUTINE  TERRKT)  and  reading  the  nine  records 
into  memory  (SUBROUTINE  INTER). 
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THE  IMAGE  GENERATING  SYSTEM 


5 , 1 Introduction 

The  components  of  the  image  generating  system  h.J  been  listed  in  Section 
1.6,  but  the  list  will  be  repeated  here  for  ease  of  reierer.ee : 

a)  For  each  of  the  four  geometry  models,  the  ray- tracing  module 
which  produces  an  unshadowed  image  of  the  subject  geometry. 

L)  For  each  of  the  four  geometry  models,  the  ray-tracing  module 
wt.jch  shadows,  with  the  subject  geometry,  a previously  gener- 
ated block  format  image  on  the  user  designated  tape  unit  IMF  I . 

c)  A scene  setting  nodule  corresponding  to  the  old  CAMERA  module 
and  performing  similar  input  reading  and  processing  functions 
for  shadow  runs  as  well  as  for  Image  generating  runs. 

d)  A data  retrieval  module  whose  components  and  functions  have 
been  described  in  Section  2. 

e)  A mooule  (PROGRAM  CONV)  for  converting  a block  format  image  to 
a string-  format  image. 

f)  A module  (PROGRAM  MERGE)  for  merging  two  block  format  images. 

g)  An  executive  whose  function  is  to  read  option  cards  and,  in 
response  to  these  option  cards,  to  invoke  the  sequence  of 
selected  runs,  while  providing  through  the  data  retrieval 
module,  the  necessary  input  for  each  run. 

Of  the  thirteen  components  named  above,  only  the  date  retrieval  module 
(a)  and  the  executive  (g)  are  completely  new.  The  vegetation  ray  tracing  of 
(a)  and  (b)  are  coeiplete  revisions  of  the  old  model  ae  described  in  Section  1.2. 
The  terrain  ray  tracing  of  (a)  and  (b)  contains  relatively  minor  revisions  as 
described  in  Sections  l.r>  and  4.1.  For  all  th#  modules  of  (a)  and  (h) , input 
and  output  have  been  altered  to  road  input  data  from  the  data  base  and  to 
output  the  block  format  usage  on  a user  designated  logical  unit  IBFO. 
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A basic  change  in  the  modules  of  (a)  allows  block-by-block  merging  of 


the  image  being  generated  with  a previously  generated  block  format  image  lo- 
cated on  the  user  designated  logical  unit  IBFI,  which,  if  set  to  zero,  will 
inhibit  the  merge.  Tha  camera-to-object  distances  for  the  points  of  a block 
of  the  established  image  are  used  to  initialize  the  maximum  search  distances 
for  the  corresponding  rays  of  the  current  image.  Thus,  there  is  a dual  ad- 
vantage in  this  procedure:  increased  efficiency  of  ray  tracing  as  well  as 
reduced  I/O  in  the  generation  of  a composite  image. 

All  of  the  modules  involved  in  the  reading  and  writing  of  block  format 
images  have  been  altered  to  conform  with  a change  in  format  introduced  to  re- 
duce I/O  time„  The  change  in  the  format  of  the  block  foraet  image  file  in- 
volves the  second  of  the  two  records  which  are  written  for  each  block  of  the 
image1.  Formerly  thi3  record  contained,  for  each  "non-sky"  point,  four  words: 
an  identifier  for  location  of  the  point  in  the  block,  and  the  product  of  the 
camera-to-object  distance  with  each  of  the  components  of  the  object  normal  at 
the  point  struck.  These  last  three  words  are  now  packed  into  a single  word, 
thus  substantially  reducing  the  length  of  the  file  for  a complex  image.  Changes 
were  thus  required  in  the  old  PICTUR  program  cf  the  .modules  of  (a)  and  in  the 
old  SHADOW  program  of  the  modules  of  (b) . In  addition  PROGRAM  CONV  and  PRO- 
GRAM MERGE  were  modified. 

The  scene  setting  components  of  (c)  are  derived  in  part  from  the  old 
CAMERA  module  and  in  part  from  the  old  SHADOW  program.  One  function  is  to 
read,  pr  corns  and  store  those  input  data  which  are  common  to  tha  elements  of 
a acenei  camera  location  and  direction,  focal  length,  image  dimensions,  grid 
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characteristics  or,  for  shadow  runs,  Bource  type,  source  direction  and  so  on. 

A second  function  is  to  determine  the  ground  plane  limits  for  retrieving  the 
trees  which  appear  in  the  scene.  On  the  other  hand  the  data  which  are  pecu- 
liar to  the  geometry  model  are  read  by  the  managing  routines  of  that  model: 
picture  number,  location  and  scale  factor  for  target*  or  i^tr.  ’•■ister  limits 
for  image  generation  or  shadowing. 

The  remainder  of  this  section  will  be  davoted  to  a discussion  of  the  execu- 
tive (Section  5.2),  a description  of  tha  overlay  structure  (Section  5.3),  a 
description  of  new  subroutines  of  tha  system  (Section  5.4),  and  finally,  a 
user's  guida  with  input  description  and  logical  unit  assignments  (Section  5.5). 


5.2  The  Executive 


The 
Figure  5 
IDB 

IPD 


IA 


IB 


S 

i 


executive' can  best  be  described  in  terms  of  the  flow  chart  of 
1 and  a glossary  of  program  variables  which  will  be  given  here: 
the  logical  unit  (set  to  1)  on  which  the  data  base 
resides. 

the  logical  unit  (user  input)  which  is  to  be  used  to 
transfer  selected  data  between  the  IDB  and  the  data 
organizing  programs,  (PREPAR)  of  the  vegetation 
nodules.  This  variable  is  also  used  as  a flag  for 
the  scene  setting  elements:  zero  for  image  runs  and 
unity  for  shadow  runs  (ISAVE=1) . 

Hollerith  data  supplied  by  the  user  on  the  first 
option  card.  Sets  the  value  of  ISAVE  according  to 
the  following  table: 


IA 

ISAVE 

Path 

ICAM 

1 

set  scene 

IGEO 

2 

' generate  an  image 

ISHA 

3 

shadow  an  image 

ICON 

4 

convert  an  image 

MERGE 

5 

merge  two  images 

IEND 

6 

rewind  tapes  and  STOP 

Hollerith  information. 

supplied  by  the  user  on  the  first 

option  card. 

Sets  the 

value  of  JSAVE  which  selects  among 

the  four  geometry  models  for  a geometry  run  (ISAVE* 2 or  3) 

and  may  serve 

to  warn 

the  scene  setting  program  (ISAVE*1) 

that  a vegetation  run 

is  contemplated. 

IB 

JSAVE 

Model 

ICG 

1 

CG  target 

I NET 

2 

camouflage  net 

IVEG 

3 

vegetation 

ITER 

4 

terrain 
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Generic  name  for  an  input  tape,  integer  value  supplied 


IBFO 


XSV1 


ISl , IS2 


JFIRST, JLAST 


by  the  user  on  the  second  option  card.  For  an  image 
generating  run,  IBFI  may  be  zero,  or  it  may  point  to  a 
block  format  image  with  which  the  current  image  is  to 
be  merged  as  it  is  generated.  For  a shadow  run  IBFI 
desingates  the  unit  of  the  image  to  be  shadowed.  For 
a CONVERT  run,  IB*  I designates  the  unit  of  the  image 
to  be  converted.  Read  as  IF(1). 

Generic  name  for  the  output  tape  unit,  user  supplied, 
for  any  run»  must  be  non-zero  for  ISAVE-2,3,4,  or  5. 

Read  as  IF (2),  except  for  ISAVE-5,  where  it  is  read  as 
IF (3) . 

Generic  name  for  a "save"  file.  If  ISV1  ? 0 , the  output 
on  IBFO  is  copied  to  I3V1,  file  number  NFILE  and  a 
printout  identifie  the  image,  the  unit  and  the  file. 
-Read  as  IF(6),  if  ISV.l  is  different  from  previous  value, 
NFILE  is  set  to  zero  and  advanced  by  unity  just  prior 
to  copy.  For  IA«IEND (ISAVE-6) , IF (I)  gives  the  sequence 
of  save  files  for  rewinding. 

Generic  names  for  units  of  two  block  format  images  to  be 

« 

merged  (ISAVE-5).  Read  as  IF(1),  IF(2).  ISl  may  also  be 
used  to  identify  a target  or  net  by  sequence  number 
(ISAVK-2,3;  JSAVT»1 , 2 ) . Then  ISl  is  lead  as  IF<4). 

The  first  and  last  indicss  of  tree  models  to  be  processed 
for  a vegetation  run,  image  or  shadow.  Read  as  IF(4)  and 
IF(0).  Default  values  (if  either  is  zero)  are  1 and  12 
respectively.  The  corresponding  treo  identifiers  are 
given  in  Section  2.1. 
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FIGURE  5.1.  FLOW  CHART  OF  TH3  SYSTEM  EXECUTIVE 
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Thus,  on  the  .first  option  card,  the  user  selects  a given  type  of  run 
{ ISAVE)  and,  for  geometry  runs  or  a scene  setting  run,  a specific  geometry 
(JSAVE).  He  must  also  select  a value  of  IPD  on  this  card  for  two  cases: 

a)  A scene-setting  run  (ISAVE)  for  images  requires  IPDOj 
for  shadows 

b)  A vegetation  image  or  shadow  run  must  have  a non-zero  value 
of  IPD. 

On  the  second  option  card,  required  for  ISAVE>1,  the  user  selects  logical 
units  for  input  and  output  and  (possibly)  for  a permanent  storage  file.  For  a 
geometry  run  for  CG  targets  or  nets,  he  also  specifies  a sequence  number  for 
the  particular  target  or  net.  For  a geometry  run  for  vegetation  he  may  specify 
a first  and  last  index  for  tree  models  to  be  processed. 

Vegetation  runs,  image  or  shadow,  are  handled  differently  from  ",  c.>:ur 
geometries.  Because  of  the  different  tree  models  which  say  occur  in  a given 
scene  the  vegetation  run  is  performed  in  a loop  on  tree  model  sequence  number 
(see  Section  2.1  for  correspondence  of  sequence  numbers,  tree  type  numbers  and 
model  identification  numbers).  The  user  may  leave  control  to  the  program  by 
leaving  IF (4)  and  IF (5)  blank.  This  will  invoke  a complete  set  of  vegetation 
runs,  i.e.,  all  sequence  numbers  found  for  this  scene  will  be  run  for  images 
and  all  even  sequence  numbers  (low  resolution  isodel)  will  be  used  in  shadowing. 
Or,  the  user  may  control  the  run  in  one  of  two  ways:  1)  he  say,  as  mentioned 

above,  select  a first  and  last  sequence  number  for  a single  exhaustive  run,  or 
2)  he  may  opt  for  several  runs  with  selected  models.  Clearly  the  decision  will 
be  based  on  the  user's  knowledge  of  the  forest  library  and  the  scene  limits. 

The  executive  calls  eleven  primary  overlays  whose  main  programs  are  given 
in  the  next  sub-section.  The  subroutines  called  by  the  executive  and  residing 
with  it  in  the  main  overlay  are  INPIC  (for  reading  geometry  dependent  camera 
dace)  and  COPYIT  (for  saving  intermediate  output  on  tape  or  disks). 
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5.3  Description  of  the  Overlay  Structure 

The  overlay,  structure  for  the  system  is  depicted  in  the  block  diagram 
of  Figure  5.2.  Each  overlay  is  represented  by  a labeled  box  which  contains 
a list  of  subroutines  headed  by  the  main  program  of  the  overlay.  The  input 
and  output  units  are  shown  on  the  alternate  block  diagram  of  Figure  5.3. 

A few  words  should  Ire  3aid  about  the  sub-program  names  to  aid  in  identi- 
fying similar  programs  in  different  overlays.  IBM  overlay  procedure  will  not 
maintain  identical  names  in  different  overlays.  The  action  would  be  to 
move  the  first  such  module  to  the  next  lower  level  and  to  eliminate  succeed- 
ing routines  of  that  name.  Hence  the  method  adopted  here  was  to  rename  in  a 
manner  suggestive  of  the  overlay  function  while  retaining  enough  of  the  generic 
name  for  easy  identification.  Thus  the  modules  PICTUR,  WTPICT,  KPICTU,  LPICTU 
are  almost  identical  but  serve  to  generate,  respectively,  combinatorial 
geometry  images,  net  or  terrain  images,  hign  resolution  vegetation  images  and 
low  resolution  vegation  images. 

The  main  overlay  (0,0)  of  the  system  contains  the  executive,  PROGRAM 
TRAFFIC,  its  subroutines  INPIC  and  COPYJT,  and  several  other  subroutines  which 
are  used  in  several  brancher  of  the  system.  Card  input  is  read  by  TRAFFIC  arid 
INPIC. 

Several  branches  of  the  program  terminate  with  the  primary  overlay : 

(1.0)  the  scene-setting  overlay,  main  program  CAMERA.  Cards 
are  read  by  CAMERA. 

(2.0)  the  data  retrieval  overlay,  main  program  DA  TRET.  The 
data  base  is  positioned  for  later  reading  or  is  accessed 
for  selection  of  data  to  be  output  on  the  IFD. 

(3.0)  the  forest  organizing  overlay  for  images,  main  program 
PKEi’AR.  The  IPD  la  read  by  PKEPAR  and  the  forest  data 
ordered  by  image  block  in  records  of  200  trees,  output 
on  TAPE 3.  The  output  TAPE10  contains  information  to 
ordar  trass  by  distance  from  the  camera.  Scratch  tapes 
li,  12  and  13  are  used. 
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FIGURE  5*  2,  BLOCK  DIAGRAM  OF  THE  OVERLAY  STRUCTURE 


VRIGID 


' " f 

'Vs  FIGURE  5.3.  BLOCK  DIAGRAM  OF  OVERLAY  STRUCTURE  SHOWING 

4*  INPUT  AND  OUTPUT  UNITS 
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FIGUPE  5.3  Continued 
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(4.0)  the  forest  organizing  overlay  for  shadows,  main  pro- 
gram PREPR.  The  IPD,  containing  forest  d-U»  and 
local  terrain  normals,  and  tha  IBFX  (to  *)«  shadowed) 
are  read  by  PREPR  to  determine  whether  the  tree 
shadow  volume  intersects  tha  object  volume  for  a 
giver,  block.  Output  tapes  3 and  10  are  similar  to 
those  frcir.  PRF.PAR,  except  that  tha  distance  array  on 
TAPE10  is  replaced  by  distance  limit*  of  the  shadow 
volume  along  tha  ray.  Scratch  tapes  11,  12  and  13 
are  used. 

(11.0)  tne  primary  overlay  for  vegetation  shadow*,  main 
program  MADOW.  MADOX  reads  tree  model  data  from  the 
IPD,  and  calls  the  routine  VSHAD  (equivalent  to  the 
old  SHADOW  program)  to  generate  shadows  for  the  trees 
represented  on  tapes  ?■  and  10.  Only  the  low  resolu- 
tion model  ray  tracing  is  Implemented  for  shadows. 

(12.0)  the  conversion  overlay,  PROGRAM  CONV,  reads  IBFI  and 
writes  on  IBFO  tha  string  format  version  of  the  input 
block  format  image. 

(13, Ci)  the  merge  overlay,  PROGRAM  MERGE,  reads  and  merges  two 
block  format  images  on  ISl  and  162  and  writes  the  re- 
sult on  IfiEO. 

The  other  branches  of  the  program  terminate  with  one  of  two  alternative 
oecondary  overlays.  In  general  the  primary  overlay  con tains  a manager,  for 
discriminating  between  the  two  secondaries  for  satisfying  a particular  option, 
and  most  of  the  routines  Uiat  are  common  to  the  two  options  served.  Each 
secondary  overlay  contains  tha  routines  peculiar  to  it*  mission  and  the  main 
program  for  that  mission.  Tha  method  of  branching  was  chosen  to  maximize  the 
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number  of  common  subroutines  in  the  primary  overlay  and  hence  to  minimize 

* * 

repetition  of  programs.  A complete  branch  is  either  an  image  run  or  a shadow 
run  for  one  of  the  geometry  models.  Each  such  branch  reads  some  card  data 
(transaction  vector  and  scale  for  targets  and  nets,  picture  number,  raster  limits 
and  geometry  data  from  the  data  base  and  may  use  an  IBF1  as  an  input  tape.  Out- 
put is  a block  format  image  on  IBFO. 

(5.0)  the  primary  overlay  for  Combinatorial  Geometry 
images  and  shadows,  main  program  CGEXEC.  CGEXEC 
is  the  manager  which  discriminates  between  images 

(5.1)  and  shadows  (5,2)  on  the  basis  of  the  value 
of  ISAVE  (2  for  images,  3 for  shadows). 

(5.1)  the  secondary  overlay  for  CG  images,  main  program 
PICT AX. 

(5.2)  the  secondary  overlay  for  CG  shadows,  main  program 

SHADOW. 

(6.0)  the  primary  overlay  for  camouflage  net  and  terrain 

imager,  main  program  NTEXEC.  NTEXEC  is  the  manager 
which  determines  whether  to  call  for  a net  image 

(6.1)  (JSAVE»2)  or  for  terrain  image  (6,2)  (JSAVE-4). 

(6.1)  the  secondary  overlay  for  net  inmges,  main  program 
CAMNKT. 

(6.2)  the  secondary  overlay  for  a terrain  image,  main 
program  LANSCP. 

(7,0)  the  primary  overlay  for  camouflage  net  and  terrain 

shadows,  is&.in  program  NTSli/U;.  NTGHAIl  disci  is;.  > .u : <;.•! 
between  net  shadows  (10,1),  (.1GAVT>2),  mid  retrain 

i dm dows  (10,2),  (JSAVK-4) . 
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(7.1)  the  secondary  overlay  for  net  shadows,  main  program 
N5HA1J  (based  on  SHADOW)  . 

(7.2)  the  secondary  overlay  for  terrain  shadows,  main 
program  TSHAD. 

(10.0)  the  primary  overlay  for  vegetation  images,  main  program 
mamace.  MAMAGE  calls  the  appropriate  secondary  overlay 
for  a hj.gh  resolution  model  (10,1)  or  for  a low  resolu- 
tion model  (10,2),. 

(10.1)  the  secondary  overlay  for  high  resolution  tree  imagea, 
main  program  VPICT1  (old  PICTUR) . 

(10.2)  the  secondary  overlay  for  low  resolution  tree  images, 
main  program  VPICT2 . 

5.4  Description  of  New  Routines 

This  material  will  be  organized  more  or  leaa  by  overlay  and  program 
branch.  We  start  with  the  main  overlay,  referring  the  reader  to  Section  5.2 
for  the  description  of  the  executive. 

Other  new  routine#  of  the  main  overlay  aret 

SURHOUTIHE  I IIJ’J r 

This  routine  reads,  processes  and  stores  geometry-dependent  camera  data 
for  image  generation:  picture  number,  location  and  acala  of  a target,  raster 
limits  for  image  scanning. 

a)  Called  from:  TRAFFIC  (for  vegetation),  PICT  AX, 

CA.MHET,  IAKSCP. 

b)  Subroutines  called:  None 

GUO  POUT  JNJ!_  I II I ’VS 

This  routine  reads,  processes  and  n tores  geometry-dependent  data  for 
vegetation  shadow  generation:  picture  numl>er  of  shadowing  geometry,  location 
and  scale  of  uhad</«iitg  giwanetry,  tauter  limits  for  shadowing. 

a)  Called  from:  TRAWC 


b)  Subroutines  celled:  None 


SUBROUTINE  COPYIT  (TP1 ,TP2 ,ISIG) 


For  ISIG*0,  COPYIT  copies  a block  '.ormat  imaga  on  TP1  to  TP2.  For  XS?G*1, 
COPYIT  copies  a string  format  image  from  TPl  to  TP2. 

a)  Called  from:  TRAFFIC 

b)  Subroutines  called:  None 

In  OVERLAY  (1,0)  the  following  subroutines  are  new  or  revised: 

PROGRAM  CAMERA 

This  is  a revised  version  of  the  old  program,  and  is  flow  charted  in 
Figure  5.4.  A "CAMERA  run"  is  an  option  which  must  be  invoked  before  a se- 
quence of  image  runs  and  before  a sequence  of  shadow  runs.  If  a vegetation 
run  is  contemplated  in  either  sequence,  the  program  parameter  IB  must  be  set 
to  IYEG  (see  Section  5.2).  Then  JSAVE  will  be  set  to  3 by  the  executive.  JTY 
(read  as  IPD  on  the  first  option  card)  is  rero  for  image  runs  and  unity  for 
shadow  runs.  For  JSAVE-3  and  any  JTY,  camera  data  (position,  direction,  focal 
depth,  focal  plana  dimensions,  grid  specifications)  must  be  available  and  is 
read  and  processed  by  CAMERA  before  the  ground  plane  limits  for  vegetation  are 
determined  by  SUBROUTINE  COVER.  For  JSAVEyi  and  JTY-0  the  camera  data  are 
read  and  processed,  but.  nothing  further  is  done  For  JSAVEyi  and  JTY»’l,  no 
camera  data  are  needed,  but  sun  data  (sourca  type  parameter  and  direction  of 
a monodirectional  source  or  position  of  a point  source)  are  read. 

Ik)  Called  from i executive,  as  SUBROUTINE  CAMERA  for 
IBM  and  aa  OVERLAY (1,0)  for  CDC. 

b)  Subroutines  called*  ANGLES,  COVER,  SHAISAT 
SUBROUTINE  COVER 

This  routine  is  called  by  CAMERA,  if  IB-IVKG,  to  find  four  corners  of  a 
ground  plana  (sea  level)  trar>«*oid  which  is  to  ba  the  saarch  area  for  trees 
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for  the  current  Scene.  COVER  calls  INDEX  to  find  the  direction  cosines  of 
limiting  rays  and  thus  establishes  corners  as  intersections  of  these  rays  with 
the  ground.  A corner  will  be  at  infinity  if  its  ray  does  not  strike  the  ground. 
If  all  four  corners  are  on  one  side  of  the  vertical  plane  through  the  camera 
facing  the  scene,  the  near  corners  are  moved  to  that  plane  along  the  sides  of 
the  trapezoid.  Some  parameters  are  stored  for  the  search  process, 

a)  Called  from!  CAMERA 

b)  Subroutines  called:  INDEX 

SUBROUTINE  SHADAT 

Called  from  CAMERA  if  JTY*1 , SHADAT  reads  sun  data. 

a)  Called  from:  CAMERA 

b)  Subroutines  called:  None 

The  routines  of  OVERLAY (2,0)  have  been  described  in  Section  2.2. 

The  routines  of  OVERLAY (3,0)  and  OVERLAY (4,0)  have  been  described  in 
Section  2.3. 

In  OVER! AY (5,0)  only  the  manager  is  new.  The  input  reading  rdutines  GENI 
and  GROX  have  been  altered  to  read  card  images  from  the  data  base. 

PROGRAM  CGEXEC 

This  is  the  main  program  of  OVERLAY (5,0) , which  contains  the  routines  common 
to  CG  image  generation  and  shadowing.  CGEXEC  receives  the  option  flag  ISAVE 
(-2  for  images,  “3  for  shadows)  through  cooaon  and  calls  the  appropriate  second- 
ary overlay. 

a)  Called  from:  executive  (as  OVERIAY(5,0) 

b)  OVERLAYS  called:  (5,1)  for  images  and  (5,2)  for  shadows 

OVERT  AY  (ri  T 1 ) contains  the  routines  specific  to  CG  imaga  generation,  main 
program  PICTAK.  PICTAK,  the  old  CG  image  manager,  has  been  altered  to  call  INPIC 
and  to  read  card  images  from  the  data  bass.  The  SUBROUTINE  PICTUR  has  been  al~ 
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tered  to  allow  merging  of  the  image  as  it  iB  generated  with  the  block  format 
image  on  a user  'designated  IbFI. 

OVERLAY (5,2)  contains  the  routines  specific  to  CG  shadowing,  main  program 
SHADOW. 

PROCRAM  S HA  HOW 

This  is  a revision  of  old  PROGRAM  SHADOW^  to  remove  reading  of  universal 
sun  and  camera  data  (now  done  by  SHADAT  and  CAMERA)  and  to  shadow  the  image  on 
IBFI  and  output,  the  shadowed  image  on  IBFO.  This  SHADOW  reads  geometry  de- 
pendent card  input:  picture  number,  origin  and  scale  of  geometry,  and  raster 
limits  for  shadowing, 

a)  Called  from:  OVERLAYS, 0) 

b)  Subroutines  called:  SETGM,  SKATES,  REED,  WRIT,  TESTI 

In  OVERLAY (6,0) , which  contains  routines  common  to  camouflage  net  and  ter- 
rain image  generation,  only  the  manager  is  new.  The  names  of  routines  used  in 
other  parts  of  the  system  have  been  changed  to  satisfy  the  IBM  requirement  of 
unique  n wing.  The  old  name  is  preceded  by  "NT"  and  the  whole  truncated  to  six 
letters.  (PICTUR  replaced  by  NTPICT,  etc.). 

P roc;  RAM  HTKXKC 

This  is  the  main  program  of  OVERLAY (6,0) . It  receives  the  flag,  JSAVE, 
through  common,  and  decides  between  camouflage  nets  (OVERLAY (6, 1) , JSAVE- 2)  and 
terrain  (OVERLAY (6, 2) , JSAVE-4). 

a)  Called  from:  executive 

b)  Overlays  called:  (6,1)  for  camouflage  nets,  (6,2)  for  terrain. 

OVERLAY (6,1)  contains  the  routinee  specific  to  camouflage  net  image  gener- 
ation, main  program  CAMNET.  The  old  CAMNET  has  been  altered  to  call  INPIC  and 
to  accept  the  input  for  restart  option  niruber  5 from  the  data  base. 
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SUBROUTINE  INSTAL 


This  routifte  reads  the  previously  processed  net  data  from  the  positioned 
data  base  (instead  of  the  old  restart  tape)  and  then  rewinds  the  data  base  file. 

a)  Called  from:  CAMNET 

b)  Subroutines  called;  None 

OVERLAY (6,2)  contains  the  routines  specific  to  terrain  image  generation, 
main  program  LANSCP.  All  routines  are  new. 

PROGRAM  IANS CP 

This  is  the  manager  program  for  terrain  images  which  calls  the  various 

« 

subroutines  in  sequence. 

a)  Called  fromi  NTEXEC  (as  OVERLAY (6, 2)) 

b)  Subroutines  called:  INPIC,  INTER,  NTPICT. 

SUBROUTINE  INTER 

INTER  reads  the  nine  records  of  processed  terrain  data  from  the  data  base 
and  rewinds  the  data  base. 

a)  Called- from:  LANSCP 

b)  Subroutines  called:  None 

SUBROUTINE  TFIND 

This  is  the  terrain  version  of  FIND  which  is  described  in  Section  4.3. 
SUBROUTINE  LOCAL 

This  routine  has  been  described  in  Section  4.3. 

OVERLAY (7,0)  contains  all  of  the  routines  for  camouflage  net  and  terrain 
shadows.  NTS  HAD  is  the  main  program  and  is  identical  to  the  PROGRAM  shadow  de- 
scribed for  OVERLAY (5, 2).  Only  the  data  reading  routine,  SKTGM(l)  has  been 
changed. 
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SUBROUTINE  NSETGM 


This  is  the  camouflage  net  yeraion  of  SETGM.  It  reads  processed  data  from 
the  positioned  data  base  trie  and  computes  and  stores  bom  parameters  for  ray 
tracing. 

a)  Called  from*  NTSHAD 

b)  Subroutines  called:  None 

SUBROUTINE  TSETGM 

Ti.is  routine  differs  from  NSETGM  only  in  the  nature  of  the  reads  from  the 
positioned  data  base.  For  shadowing  by  terrain,  only  the  first  four  records 
of  processed  data  are  needed. 

a)  Called  from:  NTSHAD 

b)  Subroutines  called:  None 

OVERLAY  (10,0)  contains  all  of  the  routines  no —on  to  high  resolution  and 
low  resolution  vegetation  image  generation,  with  the  exception  of  the  old  SUB- 
ROUTINE PICTUR  (see  discussion  of  secondary  overlays).  The  main  program  is 
MANAGE.  Many  of  the  ray-tracing  r.sv  Lut  have  not  yet  beer,  documented. 

PROGRAM  HAMACE 

This  program  reads  the  tree  model  data  into  memory  from  the  IPD,  identifies 
the  model  aa  high  or  low  resolution  and  cells  the  proper  overlay  for  ray  tracing. 

a)  Called  from:  executive,  as  (A/ERIAT (10,0)  ' 

b)  Overlays  called:  (10,1),  (10,2) 

OVERLAY (10,1)  contains,  as  main  program,  the  revised  and  renamed  subroutine 
PICTTJR,  called  here  PROGRAM  HPICT.  The  subroutines  HRAYT2  end  HRIGID  are  respect- 
ively, the  subroutines  RAYTR2  and  RIGID  for  high  resolution  ray  tracing. 

OViUU.AYd contain:!,  as  main  program,  the  revised  and  renamed  subroutine 
PICTUR,  called  her#  PROGRAM  LPICT.  The  subroutines  LRAYT2  and  LHIGID  are,  re- 
spectively, the  subroutines  HAYTR2  and  RIGID  for  low  resolution  ray  tracing, 

OVERLAY  (1)  ,0)  contains  all  of  the  routines  foe  shadowing  by  a low  rssolu- 
tlon  model  of  a tree.  The  main  program  le  KADOW. 
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PROGRAM  MADOW 


This  program  'reads  tha  trae  modal  data  into  memory  from  the  IPO  and  calls 
the  routine,  VSHAOO,  which  governs  tha  scanning  of  the  image  being  shadowed. 

If  it  proves  desirable,  in  tha  future,  to  introduce  the  liigh  resolution  model 
for  shadowing,  this  ia  the  appropriate  routine  to  use  as  a manager  to  disting- 
uish between  high  and  low  resolution  ray  tracing. 

a)  Called  front  executive,  as  OVERLAY(11,0) 

b)  Subroutines  called i VSHAOO 

OVERLAY (12,0)  is  tha  old  PROGRAM  MERGE,  described  in  reference  1,  and  re- 
vised to  read  from  the  user  designated  files  IBFI  and  ISl  and  write  the  output 
on  a user  designated  file  IBFO. 

OVERLAY (13,0)  is  tha  old  PROGRAM  CONV,  described  in  reference  1,  and  re- 
vised to  read  and  convert  the  user  designated  file  IBFX  and  output  the  result 
on  tha  user  designated  file  IBPO. 

5.5  Input  and  Logical  Unit  Assignments  for  the 
Image  Generating  System 

The  card  input  to  the  image  generating  system  is  of  two  kinds:  1)  option 

cards,  two  for  each  option  that  is  invoked,  except  for  the  camera  option,  which 
requires  only  one,  and  2)  card  input  for  scene  setting,  specific  to  each  option. 
Option  cards  are  reed  by  the  executive  and  are  used  to  select  the  program  and 
to  designate  logical  units  for  input,  output  and  permanent  storage.  The  other 
card  input  is  read  by  some  subroutine  called  by  the  selected  program. 

The  following  logical  unite  are  permanently  assigned  to  input  or  to  a 
program  and  should  not  be  used  elsewhere  except  with  greet  caution  and  knowledge 

assigned  to  the  date  besei  should  nev*»r  be  used  elsewhere. 

used  by  the  forest  date  organising  overlays,  (3,0)  and 
(4,0),  these  tepee  should  not  be  used  elsewhere  during  e 
vegetation  run. 


of  the  eyetf*rai 
TAPE! 


TAPE  3 
TAPE1 0 
TAPE 11 

TAPE 12 
TAPS 13 
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Any  other  units  named  ir.  the  program  card  of  PROGRAM  TRAFPIC  are  at  the  dis- 


posal of  the  user,  as  described  in  Section  5.2  and  in  the  input  description 
which  follows* 

OPTION  CARD  NO.  COMTEK TS  FORMAT 

1 IA,  13,  IPD  (A4,  ex,  A4,  6X,  110) 

(see  Section  5.2) 

2 IF  (X) , 1-1,6  (6110) 

(see  Section  5.2  and 

below) 

IA-ICAM 


There  are  three  possibilities* 

1)  XPD"0,  any  I?  - camera  cards  are  required. 

2)  IPO»l,  I&-IVEG  - camera  cards  and  sun  cards  are  required. 

3)  IPO-1,  IB^iVBO  - sun  cards  are  required. 


CAMERA  CARD  NO. 

CONTENTS 

FORMAT 

1 

CP,  CD,  IOPT,  NO 

(6P10, 5X,2I5) 

* CP (3)  - camera  position 

CD (3)  - camera  direction 
I0PT-0) 
aiming  point 
(IOPT-1) 

NO  - initial  ray 

spacing  in  grid 
points  (dsfsult  8) 

2 PD,  LEPP,  FPOFF  (5P10.0) 

PD  - focal  length 

LEPP  (2)  - vertical  and 
horiaontal  di- 
mensions of  focal 
plane 

3 MAX  (2Z5) 

HAX  ( 2 ) ■ maximum  numbor 
of  : n..t  itr  points 
vertii  ally  and 
horizontally 


SUN  CARD  NO 


1 IFLA,  DS,  WSS  (5X,I5,4F10.0) 

IFLA  m source  type  flag, 

■ 0,  plane  source 
- 1,  point  source 

D8  - clearance  distance  from 
imcge  (for  IDLA-O) 

WSS (3)  “ direction  from  image  to 

source  (IFLA*«0)  or  location 
of  point  source,  IFLA=1. 


IA*IGEO 


a)  Any  geometry  (any  IB)  will  require  the  following i 


IF  (1)  “ 0 

- IBFI 
IF (2)  - IBFO 
IF (6)  - 0 

« ISV1 


no  input  image  for  merging 
logical  unit  of  input  image 
logical  unit  of  output  image 
don't  save  intermediate  output 
copy  IBFO  to  ZSV1 


b)  IB  - ICG  or  INET 

IF  (4)  - IS1 

c)  IB  - IVEG 

IF (3)  • ZBF1 


IF  (4)  , IF  (5) 


ordinal  number  of  target  or  net 
in  target  or  net  file. 

if  IBFI  was  0 for  first  tree  model 
image,  then  IF  (3)  is  ths  output  file 
for  the  cumulative  image  with  the 
second  tree  model.  IBFI  and  IBFO  will 
alternate  under  program  control. 
JF1R5T,  JLAST,  first  and  last  tree 
models  to  oe  run,  aaa  Section  5.2. 


IPD 


must  be  specified. 


d)  All  image  generating  runs  require  the  following  cards: 


Image  card  no. 


contents 


PC. “-MAT 


IPIX  ■ picture  number 


CEN,SCP 


(4P10.0) 


IA-ISHA 


CEN(3)  “ locati-^  of  geometry 

origin  in  the  real  world 

SCF  - scale  factor 

blank  card 

ISTFT,  IEND,  JSTART,  JEND 

vertical  and  horizontal  raster 
limits  for  this  image. 


a)  All  shadow  runs  (any  IB)  require  (see  IA-IQSO,  above) 
IF (1 ) - IBFI 
IF (2)  - IBPO 
IF (6)  - 0 or  ISV1 

h)  IB  «*  ICG  cr  I NET 


IF (4)  - ISl 


c)  IB  - IVEG 


IP (4) , IF (5)  ’ 
IPD 


JFIRST,  JLAST 
Bust  be  specified. 


d)  All  shadow  runs  (any  IB)  require  the  sane  cards  as  image  rune,  except 
for  the  blank  card.  It  should  be  noted,  however,  that  the  raster  linits  here 
should  cover  trie  expected  shadow  area  on  the  image  being  shadowed  and  will,  in 
general,  be  different  frees  the  linite  for  the  image  generating  run.  IPIX  refers 


to  the  shadowing  geometry  und  ie  not  used  except  for  identification  on  printout. 


a 


XA'ICONV 


a)  IB  aid  IFD  are  irrelevant  to  conversion  runs. 


b)  Conversion  runs  require 


IF (1)  - IBFI 


IF  (2)  - IBFO 


c)  The  following  additional  cards  are  required: 


CONVERSION  CARD  NO. 


CONTENTS 


FORMAT 


NPICT,  ICV,  ALF,  WS (3) 


(2110,  4F10.0) 


NPICT  - highest  value  of  IPIX 

ICV  - flag:  <0  printer 

display  only 

■ 0 display  + 
tape  output 

>0  tape  output 

ALF  • fraction  of  light 

from  sun;  remainder, 
sky  light 

WS(3)  “ vector  pointing  toward 


Card  1 is  followed  by  NVICT  groups  of  cards: 


CARD  NUMBER 


CONTENTS 


FORMAT 


Nk£  « number  of  regions  for 
this  picture  number 
(**0,  for  absent  picture) 

NCOLOR(I),  1-1,  NRS 

user  supplied  color  number 
versus  region  for  this  picture 
number. 


(1015) 


Card  2 supplied  only  for  non-sero  NRS. 
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I KG 


a)  IB  and  IPD  are  irrelevant  to  merging 

b)  Merging  requires t 


IF  ( 1 ) - IS1 


logical  unit  of  first  input  image 


IF  (2)  - Ia2 


logical  unit  of  second  input  inage 


IF (3)  - IBFO 


logical  unit  of  output  image 


c)  No  additional  cards. 


IA-IEND 


a)  (IF  (I ) , Ir,l # IMAX)  are  logical  units  being  used  to  save  output.  The 


named  files  are  rewound. 
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DEMONSTRATION  GIT  THE  SYSTEM 


§ 


> 


i 


i. 

t 


All  of  the  major  opt Ion j cf  the  system  have  been  exercised  successfully 
in  test  problems  devised  by  MAGI.  In  addition,  a large  scale  demonstration 
of  the  terrain  and  vegetation  capabilities  of  the  system,  as  specified  by  the 
government,  has  been  executed.  The  resulting  photograph  is  displayed  in 
Figure  6.1. 

The  production  steps  leading  to  this  photograph  include  the  establish- 
ment and  incorporation  into  the  data  base  of  the  terrain  library  and  the 
forest  library  for  a task  site  described  by  the  government.  These  steps  will 
be  discussed  in  Section  6.1  an<l  6.2,  respectively.  The  data  for  image  genera- 
tion, shadowing  and  color  determination  will  be  discussed  in  Section  6.3. 

The  photograph  of  Figure  6.1,  at  a resolution  of  650  x 1300  raster  points, 
required  three  hours  of  computing  time  on  the  IBM  360/65  - well  within  the  pre- 
dicted times.  The  corresponding  time  for  the  CDC-66CG  is  about  one  hour.  At 
this  resolution,  and  seen  in  full  color  (which  unfortunately  is  not  possible 
in  this  report),  the  photograph  is  of  surprisingly  good  quality  and  displays  a 
realism  beyond  our  expectations. 

It  must  be  concluded  that  the  code  system  described  in  this  report  con- 
tains all  of  the  technology  required  to  produce,  with  reasonable  manpower, 
and  in  reasonable  computing  times,  high  quality  realistic  photographs  of  large 
scale  terrestrial  scenes.  The  degree  to  which  the  data  base  and  resulting 
photographe  will  represent  an  actual  test  site  will  be  limited  largely  by 
the  input  description  of  that  site. 

6.1  Terrain  Doecriptlor. 

A square  area,  1400  kilometers  on  a aide,  waa  selected  by  t:h«  government, 
and  a terrain  description  was  provided  in  the  following  form 
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a)  Elevation  values,  in  feat,  from  an  elevation  map,  for 
92  points.  A few  of  th&se  are  plotted  in  Figure  6.2, 
for  reference.  The  rest  are  given  in  Table  6.1. 

b)  A few  country  roads,  produced  from  memory.  These  are 
also  plotted  in  Figure  6.2. 

c)  Several  ground-to-ground  photographs  of  the  area,  which 
were  to  provide  some  feeling  for  color  boundaries  and 

colors . 

Ideally,  the  road  descriptions  and  color  information  should  come  from  a 
single  overhead  photograph,  with  known  camera  position  and  with  a couple  of 
landmarks  to  tie  the  photograph  to  the  elevation  map. 

Since  such  a photograph  was  not  available,  color  boundaries  were  selected 
by  MAGI  to  provide  a mix  of  sandy  and  grassy  areas  similar  to  that  of  the 
available  ground-to-ground  photographs.  These  boundaries  are  indicated  in 
Figure  6.2  by  dashed  lines.  It  was  intended  that  these  areas  be  more  sandy 
than  the  rest  of  the  regional  hence  several  of  them  were  drawn  from  points 
of  high  elevation  down  to  points  of  low  elevation,  i.e.,  to  behave  like  "wash- 
outs'’. One  of  these  regions,  namely  that  marked  12,  was  represented  as  a 
wide  road,  rather  than  as  U)  explicit  color  boundary . 

These  data  pressed  the  limits  of  the  arrays  and  the  terrain  Installation 
code  and  the  terrain  ray-tracing  codee.  fiowever,  the  computer  time  consumed 
by  the  installation  (Program  IM3TAL)  was  less  than  three  minutes  on  the  IBM 
360/65,  corresj*>nding  to  lees  than  a minute  <>n  the  CDC  6600. 

6.2  Treo_ 

The  ixmrvlar i ee  of  forested  area*,  several  lines  of  trees  and  several  in- 
dividual tm's  witfn  i ( led  l>y  the  government,  nti  (sttf>wn  in  Figure  6.3  and 

Table  (>.2.  Two  tree  typ«H  were  designated  the  pine  and  the  oak  tree  - and 

the  required  means,  variances  and  densitlas  wars  provided  by  the  government 
as  described  in  Section  J of  this  report. 
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Table  6 . 1 Terrain  Elevation  Table 
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0 
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Table  6.2  Tree  Data  For  Figujce  6.3 


REGION 

REGION  TYPE 

TREE  TYPE 

NUMBER 

HEIGHT  RANGE 

(jaeters) 

1 

Individual 

21 

1 

17. 

1 

2 

Individual 

21 

1 

14. 

8 

3 

Individual 

21 

1 

13. 

1 

4 

Individual 

21 

1 

16. 

4 

A 

Boundary 

11 

760 

2.5  - 

7.5 

21 

210 

9.6  - 

14.4 

A 

Line 

11 

113 

2.8  - 

5.2 

A 

Line 

11 

32 

0.5  - 

1.5 

B 

Boundary 

11 

82 

0.5  - 

1.5 

r 

Boundary 

11 

21 

2.0  - 

6.0 

21 

108 

2.5  - 

7.5 

E 

Line 

11 

205 

6.3  - 

7.7 

I 

Boundary 

11 

24 

1 

• 

f~4 

1.7 

'J 

Boundary 

11 

103 

6.3  - 

7.7 

II 

Boundary 

11 

471 

7.2  - 

8.8 

1 

Lina 

11 

42 

7.2  - 

8.8 

J 

Line 

11 

13 

6.3  - 

7.7 

< 

Line 

11 

40 

1.2  - 

4.5 

21 

35 

1.6  - 

6.0 

it 

Boundary 

11 

, 39 

1.0  - 

1.9 

21 

39 

1.5  - 
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M 

Boundary 

11 

17 

1.6  - 

2.4 

21 

12 

2.0  * 

3.0 

N 

Line 

11 

29 

6.4  - 

9.6 

0 

Boundary 

11 

14 

7.0  - 

13.0 

V 

boundary 

11 

40 

4.0  - 

6.0 

0 

ijoundary 

11 

52 

2.8  - 

5.2 
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Again,  an  overhead  photograph  pould  have  served  to  provide  these  inputs 
for  the  forest  generator  more  reliably  than  is  possible  from  memory. 

PROGRAM  FoRGEN  produced  something  more  than  2000  tree  boxes  in  about 
a minute  of  computing  time  on  the  IBM  360/65.  A printer  plot  (elongated)  was 
provided  by  the  program  PLOTFR  and  is  shown  in  Figure  6.4. 

6. 3 The  Photograph 

The  camera  position  and  direction  were  selected  by  the  government  and 
are  shown  in  Figure  6.1  as  CP  and  CD,  respectively.  The  position  is  near  the 
top  of  a sharp  rise  (x-450  a,  y-1200  a,  z-67.5  ■)  and  the  direction  is  slightly 
below  the  horizontal  (about  5*)  and  nearly  coincident  with  the  negative  y 
axis  (CD  - (+  0.1,  -0.99,  -0.1)).  The  sun  was  placed  behind  the  camera  with 
rays  at  an  angle  of  45*  with  the  vertical  and  along  the  y-axls . 

The  field  of  view  was  of  half-angle  22.5*  in  the  vertical  direction  and 
half-angle  45*  in  the  horizontal  direction.  Since,  the  required  resolution 
(650  x 1300  raster  points)  exceeded  the  systsm  limits  of  1000  raster  points 
in  cither  direction,  the  field  v*-?  divided  down  the  middle,  and  two  offset 
images  (650  x 650)  were  generated.  These  were  later  combined  by  a photographic 
laboratory. 

Table  6.3  is  a listing  of  the  input  required  by  the  image  generating 
eystesi  to  produce  an  image  of  the  right  hand  field  of  view.  Only  one  pair  of 
numbers  need  by  changed  (in  this  case  only  one  number)  to  produce  the  left 
hand  field  of  views  the  offset  pair  (0,-0. 5)  must  be  replaced  by  the  pair 
(0.0,  +0.5). 

The  total  running-time  for  the  shadowed  picture  has  been  mentioned  already 

as  i hr«e  hours  on  the  fOM  360/65.  About  half  of  this  was  needed  for  iaviglng 
vegetation,  somewhat  lose  than  a quarter  of  it  for  imaging  terrain  and  boqo- 


what  leae  than  a quarter  of  the  time  was  uaed  in  shadowing  with  the  vegetation. 
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6.3  Saapl*  Input  to  tim  Iaaga  Gan«r*ting  Sy«tata 
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About  ten  minutes  were  required  for  the  conversion  peas.  No  terrain  shadows 
were  possible  with  the  sun  At  a reasonable  angle,  but  a terrain  shadow  pass 
would  have  added  about  twenty  minutes  to  the  total. 

In  default  of  good  information,  colors  were  chosen  to  be  pleasing  and 

4 

"reasonable".  Use  was  made  of  the  color -mixing  technique  to  give  texture  to 
the  ground.  All  terrain  regions  and  roads  were  given  e reddish  sand  back- 
ground to  which  varying  amounts  of  green  were  edded  for  the  terrain  regions,  and 
a small  amount  of  black  was  added  for  roads.  The  machine  time  required  for  a 
color  pass  is  a few  minutes  on  the  IBM  360/65. 

Finally,  we  call  attention  to  the  fact  that  the  trees  of  region  B 
(Figure  6.3)  are  not  seen  in  the  photograph.  These  fifty  old  pine  trees 
ranged  in  height  from  0.5  meter  to  1-5  meters  end  thus,  et  the  distance  of  two 
hundred  to  four  hundred  meters,  an  individual  could  cover  no  more  than  two  or 
three  raster  points.  The  initial  coarse  grid  on  which  rays  were  fired  had  a 
spacing  of  eight  raster  points  and  it  was  tharafore  inevitable  that  most  or 
all  or  these  small,  widely  spaced,  specimens  should  be  missed.  Although,  in 
thia  case,  most  of  the  trees  ware  too  satall  to  be  seen  even  if  atruck,  a border- 
line case  can  be  treated  in  the  following  wayi  the  vertical  and  horizontal 
raster  limits  of  the  region  Involved  should  be  determined  and  tha  camera  option 
of  the  executive  should  be  invoked  with  an  initial  grid  spacing  of  4,  2 or  1 
raster  points,  followed  by  the  vegetation  image  option  with  tha  above  raster 
limits.  This  can  be  followed  by  another  "camera"  run  with  the  coarser  grid 
before  going  on  with  the  other  elements  of  the  picture. 
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APPENDIX  A - VEGETATION  ROUTINES 

INTRODUCTION 

In  addition  to  the  vegetation  ray  tracing  routines  in  the  picture  making 
system,  there  are  two  stand-alone  programs  which  are  used  to  generate  the 
tree  data  set  for  the  data  base.  The  two  Btand-alone  programs  arei 

a.  A Tree  Construction  Program  which  generates  a high  resolution 
data  set  from  card  input. 

b.  A Low  Resolution  Construction  program  which  converts  the 
high  reeolution  date  sat  into  a low  resolution  data  set. 

These  two  programs  and  the  vegetation  ray  tracing  routines  will  be  described 
here.  Thera  ere  men y routines  in  the  present  ray  tracing  routinea  and  con- 
struction programs  that  are  identical  to  those  routines  in  the  previous  multi- 
stage vegetation  model.  Details  concerning  these  routines  can  be  found  in  re- 
ferences 1 end  4. 
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THE  TF£E  CONSTRUCTION  PROT^AM 

The  construction  program  ia  written  in  six  overlays.  Its  primary  purpoat 

is  to  generate  e twc-stage  data  set  for  an  individual  vegetation  body.  Iha 

program  requires  the  same  Input  as  the  original  multi-stage  model  with  the 

* 

exception  of  one  additional  option  card  . When  used  in  the  construction  mode, 
the  program  generates  a data  file  containingi 

a.  The  location  and  dimensions  of  all  the  bodies  (cones  and  leaf 
cloud  boxes)  of  tr>e  tree  in  its  own  coordinate  system. 

b.  Three  dimensional  sort  of  these  bodies  into  a 15  x 15  x 25  cell 
structure. 

c.  The  packed  arrays  specifying  the  interna)  structure  of  the 
lead  clouda. 

d.  Miscellaneous  ray  trace  data.  These  data  era  used  by  the  high 
resolution  ray  trace  program. 


tor  dstaiis  concerning  input,  see  documentation  for  the  multi-stags 
model. 
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cutrlay  sy-iu.':-r.j»£  and  function 


OVERLAY 

ROUTINES 

FUNCTION  OF  OVERLAY 

1 

j 

(0,0) 

j 

TR  DRIVE  i 

CROSS*  1 

WRITE  2 j 

SET  CP*  ! 

i-r  2 ; 

IV* 

R hoc*.  DATA* 
tv'UIV  2 
VEGET 
ANGLES* 

CAMERA* 

TT 

REED* 

WRIT* 

BE  LONG* 

..  --  - - -j 

Overlay  is  primarily  responaible  for  calling 
other  overlays.  Contain*  routine*  used  by 
more  than  one  overlay.  Overlay  also  gener- 
ate* data  (in  subroutine  VtGET)  for  multiple 
tree*  to  b«  used  by  multi-stage  ray  tracing 
in  overlay  (5,0) 

j 

j 

j 

(0,1) 

UNPACK  ! 

REDUCE 

! 

tirxea  data  generated  by  overlay  (1,0)  to 
create  array*  containing  the  location  and 
position  of  all  the  bodies  in  tree.  Rou- 
tines u*e  the  mult.i-*tage  data  and  produce 
t)ie  required  two  stage  data  (TRC  data  and 
leaf  cloud  box  data) 

(2,0) 

SET  UR  2 
REPACK 
PACKER 
IHCLUB 

a 

BLUING  C 

Usee  data  generated  by  overlay  (1,0)  to  gen- 
erate three  dimensional  ceil  structure.  Cell 
atructvire  data  specifies  which  of  the  many 
bodies  are  contained  in  each  of  the  cells 
covering  the  tree 

skouehcing  of  kootinss  and  data  flow 

The  laboied  diagrams  illustrate  the  calling  sequence  for  overlays  <0,0), 
(1,0),  (2,0),  (4,0)  and  (S,0)  and  the  files  used  by  the  progress.  Overlay  (3,0) 
consists  only  of  routines  responsible  for  the  construction  of  the  multi-stage 
tree.  The  linkage  and  data  ara  unchanged  and  will  not  be  diagrammed.  Similar- 
ly overlay  (5,0)  contains  only  multi-stage  ray  traca  routine*  and  will  not.  be 
diag-amsed.  Further  details  concerning  the  linkage  can  be  found  in  earlier 
docuaen tat ion . 

Note  that  in  the  diagrams  setae  of  the  data  is  described  ae  data  (a) , (b) , 
(c) , or  (d).  These  refer  to  the  two  stage  data  aa  follow* « 

(a)  location  and  dimensions  of  all  bodies 

(b)  ceil  structure  arrays 

(c)  internal  struct*  re  of  leaf  cloud* 

(d)  miscellaneous  ray  trace  data 

DESCPIPTION  AND  FUNCTIONS  OF  NEW  ROOTIKgS 

Houtines  marked  with  an  asterisk  in  overlay  (0,0)  as  well  as  all  routines 
in  overlay  (3,0)  end  (5,0)  are  taken  from  the  original  multi-stage  model.  They 
are  unchanged  and  will  not  be  described.  Flow  chart's  arc  provided  for  selected 
routines. 

10r/;pAM  THRIVE 

TRTjPJVE  is  main  di lver  program.  Its  only  function*  are  reading  the  option 
card,  opening  the  files,  celling  the  starting  routines,  and  writing  out  some 

atnt let )<;,». 
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sitroutine  WRITE 2 

WRITE2  writaa  a binary  file  on  tape  1 containing  content*  of  ccscsnon  blocks, 
FG,  COLR,  KD,  KLOUC  ,»nd  RAYTR.  These  blocks  contain  ray  tracing  Information 
which  are  necessary  for  both  tha  original  multi-stage  ray  trace  and  new  two-* 

stage  high  resolution  XAYTRACE , 

SG&RiMrriNES  TT  (IS)  AND  TT2  (IN,  IM) 

These  routines  simply  compute  procassor  time  batwaen  various  program  steps 
and  print  out  the  time  with  an  indicating  message. 

VARIABLES 

IN  - Alphanumeric  massage  to  ba  printed 
IM  - Numeric  information  in  SMisaage 

snic^irriHE  f,q<itv2  (v,  h,  g,  kb,  ft) 

< 'amputee  body  data  tor  smallest  RPR  which  would  completely  contain  a given 

cone. 

VAiP.I  AIM  Jr.S 

V (1)  - cone  vertex  vector 

HO!  - oone  hoignt  vector  * 

G ( 2 ,i ) - Rip  data 
MB  - cone  bass  radius 
XT  - tone  top  radius 

;jigiM)tfTiNr:  j\  t 

Vi/.rT  l.i  ••.n»n’  tally  un<  h*n<|»«d , however,  it  <V-es  rail  v»i  hm*  'V.-rlays  for 
cell  structure  const!  ut  t ion  lased  upan  Input  option  larameters. 


si  . - ■. . i 
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PROGRAM  UNPACK 


UNPACK  is  the  main  program  cf  overlay  (1,0).  Data  describing  the  tree  is 
stored  in  coasaon  blacks  GECM2,  DTR22,  and  GEOK4  which  were  generated  by  over- 
lay (3,0).  The  data  is  stored  in  the  multi-stage  format.  UNPACK  changes  the 
data  into  a two  stage  format. 

In  the  multi-stage  format  thsre  ars  three  types  of  bodies » cores,  RPPs, 

(leaf  cloud),  and  boxes  (containing  the  entirety  of  a lower  stage).  UNPACK 
computes  the  rotation  matrix  (R) , translation  vector  (V) , and  seals  factor  (D) 
that  would  transform  body  data  for  a lower  stage  tody  into  the  proper  "reel 
world"  body  data  (complete  tree  reference  frame) . 

This  transformation  data  is  computed  from  the  stored  "copy*  data  in  the 
multi-stage  arrays.  The  resultant  two  stage  data  which  is  produced  consists 
of  the  body  data  for  each  individual  cone  or  box  (leaf  cloud  container)  ex- 
pressed in  this  single  complete  tree  reference  frame.  The  internal  structure 
of  the  leaf  cloud  RFP  in  unchanged.  The  entire  tree  is  snclossd  in  an  RPP 
with  vi'tw  vertex  located  at  the  coordinate  origin  with  the  remainder  of  the  KPS? 
situated  entirely  in  the  positive  octant. 

See  flowchart  #1 

SUBROUTINE  REDUCE  (SIZE,  NIB) 

In  l/NPATK,  data  for  each  body  is  computed  relative  to  a complete  tree  re- 
ference frame  which  is  in  reality  stage  5 at  the  multi-stage  model.  However, 
due  to  the  construction  of  boxes  over  each  4th  stage  body  (done  in  overlay  (3,0)), 
the  covering  «PP  is  larger  than  Is  necessary.  REDUCE  considers  all  of  the  boiler 

in  the  two  stage  *i«*cr  • pt  ion  and  coqputee  parameters  o£  the  mm  Host  WT  which 
would  completely  contain  nil  »f  the  bodies.  A translation  of  coordinates  J.s 
nvule  no  that,  the  RPP  .,w>v<*r tug  the  txjdle*  will  still  Maintain  the  property  of 

having  one  vertex  located  at  the  coordinate  origin. 
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VARIABLES 


SIZE  (3)  - SIZE  (1)  and  SIZE  (2)  are  x and  y dimensions  of  the  RPP  divided 
by  15.  SIZE  (3)  is  z dimension  o*  RPP  divided  by  25. 

NIR  - Total  number  of  bodies  in  two  stage  description. 

I PI  ARAM  SETUP  2 

SETUP  2 is  the  main  program  of  overlay  (2,0).  Tha  RPP  covering  the  tree 
is  subdivided  into  three  dimensional  cells  so  that  the  entire  RPP  is  covered  by 
a lb  * 15  x 25  matrix.  The  body  data  for  each  individual  cone  or  box  is  read 
from  an  intermediate  file  (TAPE  4)  one  by  one.  This  routine  (through  its  sub- 
routines) determines  which  cells  of  the  matrix  contain  any  portion  of  the  ex- 
ternal surfc.ee  of  a given  body.  It  thereby  generates  an  array  specifying 
which  bodies  are  contained  within  each  cell  of  the  covering  matrix.  Note:  if 

ary  cell  is  completely  contained  within  a given  body,  the  cell  is  deemed  not 
to  contain  the  body  (i.e.,  only  external  body  surface*  are  considered).  The 
routine  also  packs  the  cell  content  information  into  the  words  of  the  cell 
structure  array. 

Sue  flowchart  #2 


SUBROUTINE  BELONG.  C (XL,  IR,  V,  H,  R3,  RT,  CELL) 


BELONG  C computes  which  cells  of  the  covering  cell  matrix  contain  any 
portion  of  the  external  surface  of  a given  truncated  right  cone  (TRC) . First, 
the  smallest.  RFP  which  would  completely  inclose  the  TRC  ia  constructed.  From 
this  RPP,  one  determines  the  maximum  and  minimum  cell  indices  for  cells  which 
could  have  a common  volume  with  the  TRC.  Cell  boundaries  are  represented  by 
planes  of  constant  x,-  y,  or  z.  The  intersection  of  the  cones  with  each  of 
the  possible  cell  wall  planes  is  then  computed.  The  locus  of  these  intersections 
will  either  be  an  ellipse  or  some  portion  of  a conic  section  connected  by 
straight  lines  aa  shown  below* 


CEL_  WALL  P LAKES 


aide  and  Beae  or  Top 


Intersection  with  TRC 
Li do a,  Base,  end  Top 
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ini  ram  these  cell  wall  planes,  there  a re  further  intersections  with  the  two 
•Ji.iKi.nbional  cell  boundary  lattice  and  the  7SC  loci  as  shown  on  the  previous 
diagram.  All  cells  whose  boundary  contain*  ona  of  thasa  point  interaections 
datnaed  to  contain  the  TRC. 

The  routine  also  packs  the  words  for  tba  appropriate  entries  in  the  call 
structure  array  to  reflect  that  the  particular  call  oontaina  the  given  TRC 

(body  number  XR)„ 

VARIABLES 

XL  (3)  - dimensions  of  cell 
I A - body  number 
V ( 3 ) - vertex  vector  for  TRC 
H ( 3 ) - height  Vector  for  TRC 
kb  - radius  of  TRC  base 
RT  - radius  of  TRC  top 

CELL  (IS,  1 5 1 25)  - cell  structure  array 

SUBROUTINE  IHCM ID  (GM,  C,L) 

Routine  is  called  by  LEIjONG  C.  Variable*  CM  and  GL  correspond  to  values 
for  a point  intersection  between  a TRC  and  two  dimensional  cell  boundary  lattice 
on  a given  cell  wall  plana.  INCLUD  compute*  th*  indices  of  th*  call*  who** 
boundary  contains  the  point  intersection. 

fit  JHin  )UT  INK  br;  NT_  j.HTf  HU , (Ti,  A,  ft,  GK,  X.  I,  M,  GUj.  GT,X,  CJO,  CHI) 

Routine  is  culled  by  BELONG  C.  Its  function  is  bo  dataralna  whether  th* 
Ti«’  or  TRC  top  intersect,  with  a 'liven  cell  wail  plane.  If  either  the 

Oa.ni  or  ''  (>  loe*  (n'ermut,  l>T I MT  also  compute*  tho  endpoints  of  the  line  of 


R2 


inf  * react  ion 


VARIABLES 


IBT  - Logical  variables  which  is  true  whan  intersection  occurs. 

RB  - radius  of  TRC  top  or  base  (whichever  is  being  considered) 

CB  - equals  either  0 or  1.0  0 •+  baee  1.0  •+■  top 

A (3, 3)  - elements  of  a matrix  satisfying  £ A^X^Xj  **  R(t)  O^t^l 
which  represents  outer  surface  of 

cone  sides.  Note  t » 0 -*■  base  t * 1.0  -*■  top 
B (3)  - coefficients  describing  base  or  top  plane. 

GK,  K,  L,  M - set  of  parameters  describing  cell  wall  plane. 

K - 1,  L « 2,  M ■ 3 constant  x plane  at  x ■ GK 

K-2,  L-3.M-1  constant  y plane  at  y « GK 

K - 3,  L - 1,  X • 2 constant  z plane  et  * - GK 

(GLO,  QUO) I (GLI,  GMI)  - two  and  points  of  intersecting  line 

sonr-rxrrTWK  gjc' b -v,  hi,  h3,  »3{  cell,  mis,  ;-x,  xx) 

BELONG  B computes  which  cells  of  the  covering  cell  metrix  contain  any 
portion  of  the  external  surface  of  a given  box  (leaf  cloud  container).  Again 
the  bounder iee  of  cells  are  represented  by  planee  of  constant  x,  y,  or  z.  The 
Intersection  of  any  box  face  with  these  planes  are  line  segments.  The  routine 
computes  these  line  segments  (Including  their  end  points)  for  each  of  thu 
possible  cell  well  plenee.  Point  intersections  of  ths  remaining  call  well 
boundarieo  (the  cwo  dimensional  lattice  of  cell  boundaries  in  eeoh  cell  wnli 
plane)  with  the  box  plane  line  segstents  are  then  computed.  Any  cell  whose 
boundary  contains  these  point  intersections  are  deemed  to  contain  the  box  in 
gunttt  Ion, 
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su& -•  inr tx,  Y,  z,  ill,  JJJ,  cell,  XL) 

Routine  is  called  from  BELONG  B.  Variable*  (x,  y , x)  correspond  to 
values  for  a p\_int  intersection  between  a box  faca  and  the  two  dimensional 
cell  boundary  lattice  on  a given  cell  wall  plana.  IHCLUB  ccmputaa  tha  in- 
dices of  the  ceil  whose  boundary  contains  tha  point  intarsaction. 

S ' 'ft ROUTINE  ft ACKER  (K,  I,  J,  IH) 

If  core  tian  4 bodies  are  contained  in  any  call,  it  will  require  more 
than  one  word  to  list  tha  body  numbers.  PACKER  temporarily  atoraa  these  body 
n>i:ttbers  on  random  mats  storage  (RMS)  disk. 

VARIABLES 

K - cali  number  1jcK<5625 
I - present  contents  of  cell  array 
J - number  of  bodiea  in  cell 
IR-  body  number 

SUBROUTINE  REPACK 

For  all  ceils  containing  more  than  4 bodies,  body  number  ia  stored  on 
RMS.  After  all  bodies  have  bean  considered,  subroutine  REPACK  recalls  tha 
data  and  packs  tha  words  of  tha  coll  array  (Call  (15,15,25))  and  pack*  the 
overflow  data  In  an  addition  arrav  (Call  2). 

r Ron  RAM  liAThKT 

Routine  roads  information  from  three  files  (TARE 3,  TAPK4,  and  TAPX7). 
Theaa  fti«a  all  contain  different  portions  of  the  tctel  tree  data  sat.  LATSET 
then  writ**  the  enrlie  data  eet  on  one  file  TAPE). 
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RETIRED  INPUT 


Program  requires  the  same  input  as  original  vegetation  construction 
prog ran  with  the  exception  of  tha  first  card  (option  card) . 

Option  Card  (ISELEC (I),  1-1,3) 

Format  (3  (12,  IX)) 

Where 

ISELEC  (1)  - -1  No  cell  structured  data  sat 

1 Creata  cell  structured  data  sat 

3 No  printout  of  call  structured  data 
ISELEC  (2)  - 0 No  ray  trace 

1 Ray  trace  with  multi-stage  data  snt 
ISELEC  (3)  - 1 No  ray  trace  and  no  printout  of  multi- 

stage data  eat 

0 Printout  multi-stage  data  eat 

FORMAT  OF  P^Tft  OFT 

The  nut  pur.  f'i  i«  (TAPE3)  contains  two  files  each  of  one  record  and  all 
in  binary  format. 

File  1 

There  ia  a total  of  21631  words  on  file  1 in  the  following  arrays t 

Ha  ( 3000)  Contains  integer  body  data  for  body  number  I 

MA  (2*1-1)  Contains  body  type  (5-TRCi  7-BOX) 

MA  (2*1)  Contains  pointer  to  location  in  FPD  array  containing 

real  body  data 

MA  (3001-1)  Contains  region  number 

FPD  (7500)  Contains  real  body  data 

For  TRC  - 8 words  (3  for  vertex  vector t .3  for  >*aight  vector t 

iwsv  rartiuei  and  top  radius) 

For  Box  - 12  word*  (3  for  vertex  veotcri  and  3 for  each  of  the 

3 height  vectors 
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CELL  (IS,  IS  IV)  - Contain*  body  numbers  for  those  bodiaa  in  given  ele- 
ment of  the  three  dimensional  lattice.  Ml  words  are 
packed  integers.  This  is  celled  the  cell  structure 

array. 

If  there  are  1 or  2 bodies  in  the  given  cell,  word  is  packed  as  follows: 
216M2M  body  •)  ♦ 26*(l-t  body  «)  ♦ (#  of  bodies  in  cell) 


If  there  are  greater  than  2 bodies  in  given  cell,  word  ie  similarly  picked 
except  that  the  second  body  # is  replaced  by  pointer  to  locetion  in  CELL2 
array  where  remainder  of  body  numbers  can  be  found. 

CEU.2  ( 3000)  - Contains  packed  integers  with  body  numbers  for  tloae 
cells  containing  acre  than  two  bodies.  There  are  3 additional  body  numbers 
in  sect  word.  Additional  body  msebera  for  the  given  cell  ere  found  as 
necessary  in  succeeding  words  of  the  array. 


Integer  word  is 


£ (Body  •)  1024 
1-1 


XI.EAF 


XL  ( 3 ) 

NLEAF 


X dimension  of  let  stage  leaf  cloud  RPP 

Location  of  body  data  for  let  atage  oone 
in  RPD  array 

Dimensions  of  • unit  cell 

Locetion  of  RPP  data  for  let  stage  leaf  cloud 
in  fPD  array. 


File  2 


This  file  la*  1671  words  and  contain*  the  contents  of  the  following 
common  blocks: 

/V'.,'  I /<  ' I.l-/  , ftuf  I /TLiUin/  l and  / i’AYTM./ 

'Dies*  t;<  utrgm  bi.K.ke  contain  mi scallaneoce  data  used  by  • he  ray  trace  routines, 
The  structure  arid  contents  of  thaae  ctm won  blocks  are  identical  to  those  in 
the  original  nolti-stage 


OVERLAY  (2,0) 


FLOWCHART  *|  - UNPACK 


Compute  transformation*  to  takt  4th  stag*  bodies  to  comple* 
to  tree  reference  frame 


Loop  on  4th  stage  bodies 


/ It  thisX 
^4tb  stage 
body  real  7 

Vu.  coHeL 


Transform  body  data 
to  complete  tree 
reference  frame 


Oo  to  next  lower  stage 


Compute  transformation 
from  next  lower  stage  to 

complete  tree  reference 
frame 


4th  stage 
'S*.  bodies  1 


COPY 


S Lost  X* 
body  h the' 
V tooge  ? > 


. / '•  \ 

/’the  next  > 

body  a cone, 

\RPf?  or  a . 


Go  to  next 
higher  stog< 


this  the  4th 
. stage  ? / 


RPP  LEAF  CLOIJD 


Transform  leaf  clo-  i RPP 
data  to  box  data  in  tree 
reference  frame 


Transform  cone  data  in 
tree  reference  frame 


W 


FLOWCHART  2 - SETUP2 

START 

1 

Read  from  TAPE4  the  ced  * true  turn  size  and 
number  of  bodies 


Loop  over  each  bod 


Read  body  type 


ead  Box  Data 


:all  repack 


CALL  BELONGB 


Write  out  cob  struct- 
ure  on  TAPE  4 


Pack  body  number  Into  appr- 
opriate wards  of  cell  struct- 
ure  array.  If  ced  structure 

word  contains  more  than 
4 bodies  — 

CALL  PACKER 
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LOW  RESOLUTION  CONSTRUCTION  PROGRAM 

The  low  resolution  construction  program  is  a stand-alone  program  which 
converts  a high  resolution  data  sat  into  a low  resolution  data  sat.  This 
data  set  has  the  following  information t 

a.  The  location  and  dimensions  of  all  large  TRCs  making  up  the  trunk 
and  major  branch  structure  of  tha  tree. 

b.  Three  dimensional  cell  structure  description  (consists  of  sort  of 
large  TRCs  and  mean  free  path  information  for  tha  15  x 15  x 25  cell  structure) . 

c.  Miscellaneous  ray  tracing  data 

The  program  reads  the  high  resolution  data  sat  and  determines  which  of  the 
TRCs  are  to  be  retained  as  real  bodies  (basically  all  TRCs  are  retained  except 
the  1st  stage  TRCs  from  the  old  multi-stage  modal  corresponding  to  secondary 
branches) . 

The  program  reduces  tha  size  of  tha  body  data  arrays  since  box  data  for 
leaf  cloud  containers  and  1st  stage  TRCs  are  no  longer  stored. 

A new  cell  structure  is  computed.  This  array  contains  body  numbers  for 
all  of  the  remaining  TRCs.  In  addition,  this  array  contains  mean  free  path  in- 
formation for  rays  traversing  tha  calls.  This  information  is  computsd  by  con- 
sidering how  many  TRCs  or  lsaf  cloud  boxss  nave  common  volumes  with  the  cell 
in  question. 
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VEGETATION  RAY  TRACING 
High  Resolution 

The  high  resolution  ray  tracing  produces  a block  format  image  tape  iden- 
tical to  the  image  tape  produced  by  the  multi-stage  model.  The  following 
routines  are  responsible  for  the  high  resolution  ray  tracings 


PICTUR** 

RAYTR2 

MURX* 

HDEX** 

RIGID 

PINE* 

RAYTRX 

NEWGG 

CLOUD* 

RAYTRI 

FOG* 

The  routines  marked  with  an  asterisk  are  unchanged  from  the  original  mu)  ti- 
stage  model.  Those  marked  with  a double  asterisk  are  common  to  other  types  of 
ray  tracing  (i.e.,  conventional  geometry , terrain  geometry,  or  camouflage  net 
geometry).  These  routines  will  not  be  described. 

SUBROUTINE  RAYTRX  (CP,  CD,  NC) 

This  routine  is  responsible  for  assuring  that  the  proper  trees  are  in  memory 
for  the  ray  tracing.  Due  to  the  largo  number  of  trees  possible  in  the  scene,  it 
was  necessary  to  have  pre-sorted  which  trees  could  be  struck  by  rays  from  a 
given  block  on  the  block  format  image  U.pe.  This  sorted  information  is  on  e 
disk  file  (TAPE  10).  Zn  addition  geometry  information  (box  data)  for  only  a 
limited  number  of  trees  can  be  in  memory  at  any  given  time.  Hence,  the  box  data 
for  the  treee  must  be  reed  into  memory  in  groups  of  200.  This  box  data  is 
stored  on  disk  file  TAPE  3.  RAYTRX  reads  the  file  TAPE  10  and  determinas  to 
which  of  the  200  resident  trees  ray  tracing  should  be  attempted  on  the  basis  of 
the  block  indices  for  the  current  ray.  RAYTRX  will  also  read  in  anothar  group 
of  200  treee  from  TAPE  3 aa  required. 

Called  fromi  PICTUR 

Variables t CP  (3)  Camera  position  coordinates 
CD  (3)  Ray  direction 
NC  Region  strike  number 


Routines  cel led i RAYTRI 


SUBROUTINE  RAYTRl  (CP,  CD,  NO 

RAYTRl  is  responsible  for  ray  tracing  to  tha  boxes  containing  trees.  If 
a tree  box  is  struck,  the  ray  is  transformed  into  the  reference  frame  of  the 
tree  and  internal  ray  tracing  is  p«r formed  by  RAYTR2,  If  an  internal  ray 
strike  is  made,  RAYTRl  transforms  the  normals  back  into  the  reel  world  re- 
ference frame. 

Called  front  RAYTRX 

Variables!  CP  (3)  Camera  position  coordinates 
CD  (3)  Ray  direction 
NC  Region  strike  number 

Routines  called t RAYTR2,  NEWGG 

SUBROUTINE  RAYTR2  (PIS,  CP,  CD,  NC.  IMISS) 

This  routine  is  responsible  for  ray  tracing  inside  a tree  box.  A trans- 
formed ray  is  followed  through  the  cell  structure.  Tracing  is  done  to  the  geo- 
metric solids  which  comprise  the  tree.  The  trunk  and  limb  structure  is  composed 
of  Truncated  Right  Cones  (TRC)  and  the  leaf  and  twig  structure  is  contained  in 
boxes.  If  a ray  strikes  a cone,  a normal  is  computed  and  ray  tracing  ends  for 
the  ray.  If  a leaf  cloud  box  is  struck,  the  ray  is  again  transformed  into  the 
reference  frame  of  the  leaf  cloud  containers  and  internal  leaf  cloud  tracing  is 
done  by  tha  subroutines  FOG,  CLOUD,  MURK,  end  PINE.  Due  to  the  greet  number  of 
bodies  (TRC,  and  leaf  cloud  boxes)  comprising  the  tree,  there  is  a covering 
three  dimensional  cell  matrix  over  thn  entire  tree.  The  path  of  the  ray  in- 
tersects only  a few  of  the  cells  of  the  matrix  end  treeing  ie  dona  only  for 
thoae  bodies  along  the  path  of  the  ray. 

Celled  frost  RAYTRl 

Variables t DIS  - distance  to  ray  strike  from  tree  box  surface 

CP  (3)  - coordinates  on  tree  box  surface  where  ray  enters 
CD  (3)  - rsy  direction  in  tree  box  frame 
NC  - region  strike  number 

XMX8S  - logical  variable  which  is  true  of  ray  misses 
ell  internal  tree  bodies 
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Subroutines  #calledi  RIGID,  NEWGG,  FOG 

SUBROUTINE  RIGID  (IBOX,  EXIT,  START,  CLAST,  CD,  NN) 

Subroutine  RIGID  uses  the  cell  structure  array  to  determine  which  bodies 
are  contained  within  any  given  cell  of  the  covering  three  dimensional  matrix. 
The  body  numbers  are  returned  to  RAYTR2  for  ray  tracing. 


Called  froot 

RAYTR2 

Variables  s 

IBOX  (64) 

array  containing  body  numbers  of  those 
solids  (TRCs  or  leaf  cloud  boxes)  to 
be  traced  toward 

EXIT 

logical  variable  which  is  true  if  ray 
has  exitad  the  tree  box 

START 

“ 

distance  ray  has  traveled  point  of  entry 
into  tree  box 

CLAST  (3) 

- 

coordinates  on  tree  surface  where 
ray  enters 

CD  (3) 

- 

ray  direction  in  tree  box  frame 

NN 

- 

number  of  body  nvmbers  in  array  IBOX 

Subroutines  cal led i None 

SUBROUTINE  NEWGG  (ISIDE,  PIS,  BOX,  I TYPE,  CD,  CP,  IMISS) 

This  routine  computes  the  strike  distance  for  intersections  of  rays  with 
the  three  dimensional  solids  (TRCs  or  boxes). 


Called  fromi 

RAYTRl , RAYTR2  1 

Variables i 

ISIDE 

integer  specifying  which  surface  strike  occurs 
(i.e.,  fist  surface  or  cons  side) 

DIS  . 

distance  from  ray  origin  to  striks  point 

BOX  (12) 

array  containing  body  data  (i.e.,  box  height 
vectors  or  TRC  data) 

ITYPE 
CD  (3) 
CP  (3) 
IMISS 

integer  specifying  body  type  (TRC  or  box) 
ray  direction 
coordinates  of  ray  origin 

logical  variable  which  is  true  of  ray  misses  body 

■:  ?**»&(:'  -»'-  J>v  j .-;  ■->'•■'■:  r-  ~ '- 


Routines  called;  None 
Low  Resolution 

In  the  high  resolution  model,  all  ray  strikes  are  computed  as  intersec- 
tions with  three  dimensional  solids  or  as  strikes  with  leaves  or  twigs  inside 
leaf  cloud  containers.  In  the  low  resolution  model,  ray  strikes  are  deter- 
mined probabilistically  for  leaf  clouds  and  the  smaller  cones  representing 
small  limbs.  For  the  trunk  and  major  branch  structure,  the  ray  tracing  is 
done  as  it  is  in  the  high  resolution  model.  In  the  cell  structure  array  (the 
three  dimensional  matrix  covering  the  tres),  mean  free  paths  for  rays  tra- 
versing the  cell  are  stored  instead  of  body  numbers  for  leaf  cloud  boxes  or 
small  cones.  The  body  numbers  for  the  trunk  and  major  branch  structure  are 
also  stored  in  the  cell  structure  array. 

The  routines  which  are  changed  in  the  low  resolution  model  are  RAYTR2  and 
RIGID.  In  addition,  the  routines  FOG,  MURK,  PINE,  and  CLOUD  are  no  longer  used 
since  there  is  no  ray  tracing  internal  to  the  leaf  cloud  containers. 

SUBROUTINE  RAYTR2  (PIS,  CP,  CD,  NC,  IMISS) 

The  routine  arguments  and  definition  of  these  arguments  are  the  same  as 
for  the  high  resolution  RAYTR2.  The  low  resolution  RAYTR2  trace*  only  to  those 
TRCs  comprising  the  trunk  and  major  branch  structure.  Leaf  strikes  and  small 
twig  strikes  are  computed  probabilistically.  As  the  ray  is  traversing  the 
three  dimensional  cell  structure,  the  path  length  of  the  ray  in  each  of  the 
intersected  cells  is  computed.  Stored  in  the  cell  structure  array  are  mean 
free  paths  (A)  for  rays  in  that  given  cell.  A random  number  is  drawn  (P)  be- 
tween 0.  and  then  a distance  R - - A in  (1-Pi  is  computed.  This  i*  the 
distance  a ray  would  travel  in  a medium  of  mean  free  path  A with  a probability  P 
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If  this  distance  is  less  than  the  path  length  of  the  ray  in  the  given  cell, 
then  a ray  strike  is  determined  to  be  made  at  distance  R.  Otherwise  the  ray 
is  determined  not  to  have  struck  any  leaves  and  continues  on  to  the  next  coll 
provided  no  TRC  hits  are  made.  Thera  is  never  any  ray  tracing  either  to 
leaf  cloud  boxes  or  internal  to  these  boxee. 

Called  from:  RAYTR1 

Variables:  Same  as  for  high  resolution  RAYTR2 

Routines  called:  RIGID,  NEWGG 

SUBROUTINE  RIGID  (IBOX,  EXIT,  START,  CLAST,  CD,  NN) 

This  routine  extracts  both  the  mean  free  paths  and  body  numbers  for  the 
large  TRCs  from  the  cell  structure  array  for  each  of  the  celle  along  the  path 
of  the  ray.  This  information  ie  returned  to  RAYTR2  for  processing. 

Called  from:  RAYTR2 

Variables:  - Same  as  for  high  resolution  RIGID 

Routines  called:  None 


